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Abstracts 
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A DISK MODEL FOR THE 10-MICROMETER SPEC- 
TRA OF YOUNG STELLAR OBJECTS WITH DUST- 
DOMINATED CONTINUA. J. E. Bowey 12 and A J 
Adamson 1 - 3 , University of Central Lancashire, Preston, UK, 
2 Queen Mary and Westfield College, London, UK (j.e.bowey@ 
qmw.ac.uk), 3 UK Infrared Telescope (UKIRT), Joint Astronomy 
Centre, Hilo HI, USA (a.adamson@jach.hawaii.edu). 

Introduction: Young stellar object (YSO) spectra vary 

widely [e.g., 1,2] and are characterized primarily by nonphoto- 
spheric dust emission and frequently with silicate emission or 
absorption features. The 10-pm disk model was initially developed 
to estimate the quantity of optically -thick and optically-thin emis- 
sion in the spectra of Taurus-Elias 7 (Haro6-10) and HL Tau but 
unlike previous representations, which reproduce the observed 
spectrum by modeling the flux emitted by dust as power laws, this 
model reveals information about the YSO disks themselves. The 
derived absorption profiles may be used to study the silicates and 
to diagnose the presence of ices. 

Sources: HL Tau and Elias 7 (Haro6-10) are two bright pre- 
main-sequence T Tauri stars in the Taurus Molecular Cloud 
(TMC) whose lO-pm spectra exhibit negligible photospheric 
emission, no PAH emission, and deep silicate absorption features 
(Fig. 1). The HL Tau system has been well studied at a variety of 
wavelengths [4]; its photosphere is obscured by dust associated 
with expanding cavities perpendicular to a stable 150-AU disk 
inclined at 67°. Elias 7 is a binary system with optical and infra- 
red components [e.g., 5,6] that may be surrounded by a disklike 
circumbinary envelope [5]. 



w 


Elias 7 (Haro6-10) HL Tau 



Fig. 1. 10-jim flux spectra of Elias 7 and HL Tau with modeled fluxes 
(solid) and derived continua (dashed). The derived continua of HL Tau 
are partially optically thin, with silicate emission features; those of Elias 
7 are optically thick. (Since the southern component of Elias 7 is not 
signifi -cantly extincted by silicates, the flux from this component was 
subtracted from the observed spectrum before fitting the silicate feature.) 
8-13-pm spectroscopy was carried out using CGS3 at UKIRT. 


Model: After Adams et al. [7] and Mannings et al. [8], the 

dusty disk is assumed to have radial temperature and density dis- 
tributions 


T(r) = T*(R./r)q (1) 

I(r) = £*(R*/r)P (2) 


where T* is the temperature of dust at the stellar radius, R. and 
the surface density of material, X(r) is equal to I at the stellar 
radius, q is the index of the temperature distribution, and the mass 
density index, p = 7/4, is adopted from protostellar theory [e.g., 7]. 

To avoid making assumptions about R D and unlike [7] and [8] 
we integrate over temperature instead of radius. The flux emitted 
by silicates in the disk becomes 

Fvdisk = af° B V (T)[1 - exp(-ce v r(T)~ p )]r(T) q+2 dT (3) 

J *sil 

where r(T) is obtained from equation (1). 

Assuming that optically thin emission in the 10-pm band is 
dominated by silicates, the emission optical depth, x v (r,0), is given 
by 


x v =ks v 


So 


cosO 


^T=ce v rOT 

r 


(4) 


where 0 is the angle of inclination of the disk to the sky plane and 


c k£ 0 Rg 

COS0 


(5) 


k is a constant that scales the silicate emissivity, £ v . 1 and R 0 are 
defined as the mass density and radius where T(r) = T 0 = 1000K 
at the inner edge of the silicate disk. is estimated from equa- 
tion (1), with R„ = 2R e , the radius of a typical T Tauri star [8]. 

The intrinsic flux of the disk is extincted by material in the 
cooler outer regions of the disk and foreground molecular cloud 
the flux observed, F v , becomes 

F v = F vdisk exp{-t cont (X) }exp(-c sil e v } (6) 

where e v (normalized to unity at its peak) is assumed to be the 
same within the protostellar disk and TMC. The continuum opti- 
cal depth, T cont (X) is derived from A v using a power-law extinc- 
tion curve based on that observed in the p Oph Molecular Cloud 
[see 3]. 

Parameters a, q, c, and c sil are fitted by ^-minimization. The 
silicate emissivities used during fitting were those that best repre- 
sent the silicate absorption feature of dust in quiescent regions of 
the TMC toward the field star Taurus-Elias 16 [3]. The quality of 
fit is insensitive to changes in T* and A v . However, changes in T* 
and A v affect the derived disk radius. 

Fits and Derived Absorption Profiles: The flux spectra of 

Elias 7 and HL Tau, like those of the TMC sources Elias 16 and 
18 [3], are best fitted with the Trapezium emissivity [9] {X\ = 24 
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and 9. 1 respectively); the flux spectra, fits and derived continua are 
given in Fig. 1. Structure in the modeled continuum of HL Tau 
indicates that the disk emission is optically thin in regions of the 
disk that emit significantly at 10 pm; Elias 7 is optically thick. 

The derived absorption profiles of HL Tau, Elias 7 and Elias 
16 (a field star obscured only by cold molecular cloud material 
[3]) are shown in Fig. 2; the similarity of these profiles [10] sug- 
gests that the 10-pm disk model has successfully separated the 
10-pm emission and absorption components. 

Radial Temperature Distribution: The fitted temperature in- 

dices, q = 0.33 and 0.43. for Elias 7 and HL Tau, respectively, are 
lower than the 0.5-0.61 range obtained by [8] using the ALS disk 



8 9 10 11, 12 13 
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Fig. 2. Normalized 10-pm optical depth profiles of TMC dust: Elias 16 
[3], quiescent molecular cloud; Elias 7 and HL Tau, YSO and foreground 
molecular cloud. Each spectrum is offset by 0.5 from the one below it. 


model [7]. However, our results are consistent with the MIR tem- 
perature distribution produced by Boss’ [11] radiative and hydro- 
dynamical protoplanetary disk models for T Tau [12]. In these 
models the radial power-law temperature distribution is modified 
in the disk mid-plane (1 AU < r < 10 AU), by the effects of grain 
condensation and evaporation on opacity, causing T(r) to plateau 
at radii of -2-4 AU, which corresponds to the 200-300 K tem- 
perature range. 

Inferred 10-pm Disk Radius, R sa : We may infer the radius 
of a 10-pm image using equation (1). We define the radius of the 
disk when k = 10 pm, R sil , to be that where the interpolated an- 
nular flux, F^^ = F pk /e. F^ is the largest annular flux compo- 
nent at 10.0 pm. The -400 AU 10-pm disk of Elias 7 (T* = 4000 
K; A v 6-10 m ) is consistent with the extent of the flattened 
circumbinary envelope proposed by [5], In the L band the radius 
of our disk would be about 20 AU, or -0."14, which roughly co- 
incides with the 75% contour of their L-band residual emission. 
However, our estimate of 730-1400 AU (T* = 4400 K; A v = 12- 
22 m ) at 34 K for HL Tau is not consistent with the value of 
100 AU obtained by [13] using the ALS model. Given the non- 
uniformity of the mass-density and temperature distributions with 
radius, reliable estimates of the radius and mass of the entire 
infrared-submillimeter disk cannot be obtained from fits with the 
10-pm disk model. 

References: [1] Hanner M. S. et al. (1995) Astrophys . J., 

438, 250-258. [2] Hanner M. S. et al. (1998) Astrophys. J., 502 , 
871. [3] Bowey J. E. et al. (1998) MNRAS, 298, 131. [4] Close 
L. M. et al. (1997) Astrophys. J., 478 , 766-777. [5] Menard F. 
et al. (1993) Astrophys. J. Lett., 414, L117-L120. [6] Van Cleve 
J. E. etal. (1994) Astrophys. Space Sci, 212, 231-238. [7] Adams 
F. C. et al. (1998) Astrophys. J., 326, 865-883. [8] Mannings V. 
and Emerson J. P. (1994) MNRAS, 267, 361-378. [9] Forrest 
W. J. etal. (1975) Astrophys. J., 195, 423-440. [10] Bowey J. E. 
(1998) Interpreting the 10-pm Astronomical Silicate Feature, 
Univ. of Central Lancashire, Preston, UK. [11] Boss A. P. (1996) 
Astrophys . J., 469, 906-920. [12] Boss A. P. and Yorke H. W. 
(1996) Astrophys. J., 469, 366-372. [13] Adams F. C. et al. 
(1990) Astrophys. J., 357, 606-620. 

THE MINERALOGICAL EVOLUTION OF SILICATE 
DUST IN INTERSTELLAR ENVIRONMENTS. J E 

Bowey 1 * 2 , A. J. Adamson 13 , and A. K. Speck 14 , University of 
Central Lancashire, Preston, UK, 2 Queen Mary and Westfield 
College, London, UK (j.e.bowey@qmw.ac.uk), 3 UK Infrared 
Telescope (UKIRT), Joint Astronomy Centre, Hilo HI, USA 
(a.adamson@jach.hawaii.edu), 4 University College London, 
London, UK (aks@star.ucl.ac.uk). 

Introduction: Attempts to fit the 10-pm astronomical silicate 

feature in a wide variety of environments (circumstellar. diffuse 
medium, YSO and comets) using the spectra of single specimens 
of amorphous laboratory silicate have been largely unsuccessful 
[e.g., 1-3]. In contrast to previous assumptions, 20-45-pm ISO 
spectra indicate the presence of crystalline silicates in circumstel- 
lar dust and the Red Rectangle [4,5] and grains from meteoritic 
sources (GEMS, crystalline IDPs, and chondrules) exhibit signifi- 
cant chemical and mineralogical variation. Similarly we find that 
the average 10-pm spectrum obtained from a number of finely 
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determined 250 K-emission profiles of O-rich circumstellar dust 
(where the UV flux is insufficient to excite PAH or ionized PAH 
bands if a carbonacous dust component were found) and in the 
absorption profiles of dust in the diffuse medium and quiescent 
and YSO environments of the Taurus molecular cloud (where the 
abundances of HMT and methanol ice are consistent with zero 
optical depth [2]) are associated with those of crystalline miner- 
als. Absence of a distinct 11.2-pm feature in interstellar spectra 
suggests that the fraction of crystalline olivine in these environ- 
ments must not be large, while the blending of bands due to mixed 
crystalline pyroxenes makes their identification in 10-pm spectra 
difficult. However, narrow structures at 9.4 pm (TMC YSOs) and 
9. 7-9. 8 pm (circumstellar dust) indicate the presence of layer 
silicates similar to talc and clay minerals, and together with a com- 
mon serpentine (10.0-10.6 pm) band suggest that hydroxylated 
minerals with spectral properties similar to terrestrial layer silicates 
are a ubiquitous component of cosmic silicate dust. 

Evolution: The silicate profile peaks at a wavelength that 

decreases between circumstellar, diffuse medium, and molecular- 
cloud environments, indicating [after 7] that the amorphous pyrox- 
ene content of initially olivine-rich interstellar dust increases 
between these environments, and hence with the age of the grains. 
In geology, direct olivine-plus-Si0 2 to pyroxene reactions occur 
only at high pressure within the terrestrial mantle and Fe-rich 
silicates are extremely resistant to chemical alteration while Mg- 
rich olivines are among the most reactive of silicate minerals. 
Therefore the most plausible mechanism for the increase in pyrox- 
ene is the hydroxylation of Mg-rich olivine to form a serpentine- 
like silicate that is subsequently annealed to become a mixture of 
amorphous pyroxene and olivine. This mechanism is particularly 
appealing, since terrestrial and laboratory olivine samples are 
readily converted to serpentine in the presence of water and, af- 
ter extended annealing of the hydrated silicate, the first crystalline 
band to appear is the 1 1.2-pm olivine feature frequently observed 
in cometary spectra. 

References: [1] Mutschke H. et al. (1998) Astron. Astrophys., 

333, 188. [2] Bowey J. E. et al. (1998) MNRAS, 298, 131- 
138. [3] Jaeger C. et al. (1994) Astron. Astrophys., 292, 641-655. 
[4] Waters L. B. F. M. et al. (1996) Astron. Astrophys 315 , 
L361-L364. [5] Waters L. B. F. M. et al. (1998) Nature , 39, 868- 
871. [6] Hanner M. S. et al. (1994) Astrophys . J., 425 274-285. 
[7] Guertler J. and Henning T. (1986) Astrophys. Space Sci., 128, 
163. [8] Ferraro J. R. (1982) The Sadder Infrared Spectra Hand- 
book of Minerals and Clays, Sadtler Research Laboratories. 

[9] Bowey J. E. and Adamson A. J. (1999) JGR, submitted. 

[10] Bowey J. E. (1998) Ph.D. thesis. Interpreting the 10-pm 
Astronomical Silicate Feature , Univ. of Central Lancashire, 
Preston, UK. [11] Bianco A. et al. (1993) ASP Conf. Series, 41, 
267. [12] Speck A. K. et al. (1999) in preparation. 


COMPARISON OF COLLECTED INTERPLANETARY 
DUST PARTICLES WITH DUST IN SPACE. J P Bradley 
and L. P. Keller, MVA Inc., Norcross GA 30093, USA (jbradley@ 
mvainc.com). 

Infrared (IR) spectroscopy enables direct comparison of dust in 
space with interplanetary dust particles (IDPs) in the laboratory 
[1,2]. The Infrared Space Observatory (ISO) has provided new' 
data about circumstellar, interstellar, and cometary dust. One of the 


most important findings of ISO is the detection of crystalline Mg- 
rich silicates (together with glassy silicates) around both young and 
old stars and Comet Hale-Bopp [3-5]. (Comet Halley’s dust also 
contains Mg-rich silicates [6].) The ISO findings are significant 
because Mg-rich silicates are rare in primitive meteorites (and 
polar micrometeorites) but abundant in IDPs, specifically the fluffy 
chondntic porous (CP) subset of IDPs, some of which are almost 
certainly from comets [7]. Together with glassy silicate grains 
(i.e., GEMS), submicrometer enstatite (MgSi0 3 ) and forsterite 
(Mg 2 Si0 4 ) are the major crystalline silicates in CP IDPs. Some of 
these crystals exhibit morphological, crystallographic, and/or com- 
positional evidence of vapor phase growth [7]. It has been pro- 
posed that GEMS are surviving presolar interstellar “amorphous 
silicates” [6]. If true, the IR silicate features should be similar to 
those observed in astronomical spectra. 

Previous studies have demonstrated that it is possible to mea- 
sure the IR spectral properties of individual IDPs and thin sections 
of IDPs [1,2]. Obtaining IR data from “pure” GEMS is exceed- 
ingly difficult because (1) they are small with subnanogram 
masses, (2) they are usually intimately mixed with other silicates 
(e.g., enstatite and forsterite), and (3) they absorb weakly in the IR 
(relative to crystalline silicates). However, the IR beamline at the 
National Synchrotron Light Source (NSLS) at Brookhaven Na- 
tional Laboratory is uniquely suited for IR spectroscopy of GEMS 
because it provides a spot >100x brighter than conventional labo- 
ratory IR benches [9], plus a detector covering the 2.5-25-pm 
(400-4000 cm -1 ) wavelength range. 

We have collected IR spectra at NSLS from GEMS-rich re- 
gions of thin sections of IDPs sectioned in a sulfur and supported 
on carbon-film TEM grids, GEMS-rich fragments of IDPs 
mounted on carbon-film TEM grids, and whole (CP) IDPs pressed 
into potassium bromide (KBr) substrates. IDPs mounted and sec- 
tioned in sulfur (rather than epoxy) are most suitable for the IR 
measurements because the sulfur can be removed (sublimed) 
after sectioning, thereby minimizing spectral background in low 
signal -to-noise GEMS spectra. 

Transmission IR spectra are usually collected over the wave- 
length range of 4000-400 cm* 1 (2.5-25 jam) with an aperture (up 
to 20 x 20 pm) positioned around each specimen. An open area 
(no sample) serves as a reference. Background spectra are acquired 
from the carbon thin-film substrate alongside the thin section. 
Final baseline corrected (and smoothed) spectra are obtained by 
subtracting the background from the sample spectrum. Between 
256 and 1000 scans (interferograms) were collected from each 
specimen with spectral resolution of 4 cm -1 . 

Most “GEMS-rich” thin sections also contain significant 
amounts of submicrometer enstatite and/or forsterite crystals. Their 
10-jim (Si-0 stretch) features are broad (2-3.5 pm FWHM) with 
absorption maxima between 9.3 and 10 pm and often show a 
broad excess absorption on the long-wavelength wing of the pro- 
files. If a thin section contains significant amounts of enstatite, the 
absorption maximum trends toward -9.3 pm. If there is significant 
forsterite, the absorption maximum shifts toward 10 pm with a 
pronounced shoulder at -11.2 |im. Preliminary spectra have also 
been collected from regions of essentially “pure” GEMS regions 
of several thin sections. In all cases, the 10-jim silicate features 
exhibit absorption maxima at -9.8 pm. 

We have also measured the 10- (Si-O stretch) and 18-pm (Si- 
O bending) features of two chondritic porous IDPs pressed into 
KBr. Both IDPs exhibit glassy “GEMS-like” silicate features with 
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Fig. I. Unnormalized spectra of Mg-rich, Fe-rich, and Si0 2 minerals — 
layer silicates (clays, talc, and serpentines; solid); pyroxenes or chain 
silicates (dashed); olivines (dotted); and silicas or Si0 2 minerals (dot- 
dash). The strongest band is that due to talc. Olivine, pyroxene, and 
serpentine spectra are due to individual members of each solid solution 
series in which the ratios of one metal to another vary. 


structured crystalline minerals of cosmic elemental composition is 
an almost smooth, broad band (Fig. 1). Fine structure in interstel- 
lar profiles is closely matched by an average spectrum obtained 
from a subset of crystalline MgFe and MgFe-rich meteoritic min- 
erals with or without a fraction of amorphous silicate. 

Observations: Molecular cloud absorption profiles are 

broader than those of circumstellar and diffuse medium silicates, 
due mainly to a ~2-pm-wide excess centered at ~1 1.5 pm (Fig. 2). 



Wavelength/ /tm 

Fig. 2. The 10-pm silicate feature in astronomical environments (solid 
with error bars); average of plotted crystalline minerals excluding silicas 
(dashed); 1:1 average of crystalline and amorphous silicate (dot-dash); 
average of five amorphous (KRS-5 substrate) spectra from [11] (dotted). 
Spectra are normalized to unity and offset in the y-axis by multiples of 
0.5. References: circumstellar spectra [12]; Elias 16 and Cyg OB2 no. 12 
[2]; TMC (Taurus molecular cloud) + YSO [9,10]. 


This shoulder has previously been interpreted as the shifted peak 
characteristic of large grain sizes (<1 pm) [e.g., 2,6], however, it 
could indicate the presence of a large fraction of chemically-var- 
ied crystalline pyroxenes or the existence of more Al-rich silicates 
in molecular clouds. 

Fine Structure: Bearing in mind the uncertainties due to 

telluric ozone, wavelength shift with laboratory matrix, gram size, 
and the metal content of crystalline silicates, fine structure in well- 
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broad 10 jam profiles (2.1 and 2.9 FWHM respectively), absorp- 
tion maxima at 9.2 pm and 9.8 pm, broad excess absorption on 
the long- wavelength wing of the profiles, and an 11.2 pm “bump” 
due to crystalline olivine. (Weak features at — 1 6.5 pm may also 
be due to olivine.) The -18-pm features are broad and featureless 
with absorption maxima at 18.9 pm and 19.4 pm respectively. 

How do astronomical silicates compare with GEMS? First, 
10- pm astronomical silicates exhibit significant variation. For ex- 
ample, there is a decrease in the wavelength of maximum absortion 
(or emission) between circumstellar, diffuse medium, and molecu- 
lar cloud environments that has been attributed to the “amorphous 
pyroxene” content of initially olivine-rich dust with age [10,11]. 
Features from individual GEMS-rich specimens exhibit similar 
variation. The similarity is intriguing because attempts to fit the 
10-pm astronomical silicate feature in a wide variety of environ- 
ments (circumstellar, diffuse medium, young stellar objects, and 
comets) using the spectra of other amorphous and crystalline sili- 
cate minerals have been largely unsatisfactory [12-14]. The great- 
est discrepancy occurs in the FWHM of the features since silicate 
minerals range between 1.7 and 2.5 pm, while the astronomical 
features are typically >2.5 pm wide. Molecular cloud profiles (e.g.. 
Trapezium) are broader than those of circumstellar and diffuse 
medium silicates, due mainly to a -2-pm-wide excess centered at 
-11.5 pm. This feature has been interpreted as the shifted peak 
characteristic of large grain sizes (>1 pm) [13,15,16], Or indica- 
tive of the presence of a large fraction of compositionally variable 
silicates [10]. The ~18-pm astronomical feature is broad and fea- 
tureless, peaks at -18.5 pm, and is relatively constant. The two 
IDP ~18-pm features so far measured peak at significantly longer 
wavelengths. 
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INTERPRETATION OF EXTRAGALACTIC EXTINCTION 
MEASUREMENTS USING THE MAXIMUM ENTROPY 
METHOD. G. C. Clayton 1 , M. J. Wolff 2 , K. D. Gordon 1 , K. A. 
Misselt 1 , department of Physics and Astronomy. Louisiana State 
University, Baton Rouge LA 70803, USA (gclayton@fenway. 
phys.lsu.edu), 2 Space Science Institute, Suite 294, 1234 Innovation 
Drive, Boulder CO 80303-7814, USA (wolff@colorado.edu). 

Introduction: Understanding interstellar dust and its role in 

the universe are very important. Virtually all observations of as- 
trophysical objects and their physical processes are affected by the 


presence of dust either within the system being studied or along 
its line of sight. Properties of galactic dust have been amenable to 
study through observations of scattering and absorption along 
individual sightlines toward point sources, and through high spa- 
tial resolution observations of extended sources. The assumption 
of galactic dust properties in modeling the effects of extragalac- 
tic dust can and does lead to large errors. With the aid of NASA 
observatories such as Hubble Space Telescope (HST), interstellar 
dust in the Local Group galaxies can be subjected to the same kind 
of observations used to study galactic dust. Individual sightlines 
can be observed at high spatial resolution allowing simple geom- 
etries to be studied. 

Results and Discussion: We present size distributions of 

interstellar dust grains derived from several sightlines in three 
Local Group galaxies: the Small and Large Magellanic Clouds, 
and M3 1 . Employing the maximum entropy method (MEM) algo- 
rithm, as developed by Kim et al. [1-3], we fit interstellar extinc- 
tion measurements which span the wavelength range 0.125-3 pm. 
In order to facilitate comparison, our sample includes several ga- 
lactic sightlines whose extinction curves share some similarities 
with the extragalactic targets. Our work has greatly benefited from 
the recent reanalysis of IUE observations of the Magellanic Clouds 
[4,5]. The grain models employed for this presentation are the sim- 
plistic homogeneous sphere models (i.e., [6] (MRN)) with two 
(graphite, silicate) or three (graphite, silicate, and amorphous C) 
components. Though such usage is only a first step, the results do 
provide interesting insight into the use of grain size as a diagnostic 
of extragalactic environment. 

The versatility of the MEM technique allows one to easily 
include a range of grain models in the analysis, from the simple 
MRN-like approach, to more realistic composite grains that can in- 
clude the effects of porosity, solid inclusions (silicate, graphite, 
superparamagnetic, etc.), and alignment. 
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BUILDING A NETWORK OF SELF-CONSISTENT MID- 
INFRARED CALIBRATION STARS. M Cohen 1 , R G 
Walker 1 , and F. C. Wittebom 2 , Vanguard Research Inc., 5321 
Scotts Valley Drive, Suite 204, Scotts Valley, CA 95066, USA, 
2 Mail Stop 245-6, NASA Ames Research Center, Moffett Field 
CA 94035-1000, USA. 

Introduction: Once spectrometers are operating under stable 

reproducible conditions (e.g., in space), the ultimate limitation to 
spectral accuracy and reliability becomes our knowledge of cali- 
bration stars. I discuss the methods we have used to establish an 
all-sky network of absolute calibrators spanning a 12-pm dynamic 
range of 160 [1]. We have taken great care to assure the self- 
consistency of this network: that is, whether using a faint or a 
bright standard, one obtains the same pV/s detector output per 
W/cm 2 incident irradiance. 

Network: The network consists fundamentally of an inverted 

pyramid. At the apex are two hot stars represented by Kurucz 
theoretical models (Sirius and Vega) [2]. The second rank consists 
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of 12 bright KO-MOIIIs; and derived from these are 422 fainter 
cool giants within the same spectral range. There are direct ties to 
the absolute calibrations of COBE/DIRBE [3], both the near- and 
mid-infrared IRTS spectrometers; MSX; ISO, particularly the 
SWS; and WIRE. First space-based validations from the IRTS/ 
NIRS demonstrate this self-consistency, while MSX absolute ra- 
diometry (analysis still in progress) assures a precision of better 
than 3% for the tier of bright giants and Sirius, our primary ref- 
erence star. 

Future Developments: All the observed spectra (as distinct 

from the Kurucz models) are low-resolution (resolving power 
-50-250), 1.2-35 pm complete (except for tiny gaps in the ter- 
restrial C0 2 bands at 4.3 pm). I describe how one might extend 
these calibrators down in intensity and up in resolving power to 
meet the needs of future groundbased, airborne, and spacebome 
spectrometers. Both these advances will be necessary to support 
the new generation of infrared instruments and astronomical ques- 
tions. 

Goals: One long-term goal is to use theoretical model atmos- 

pheres to predict the emergent spectra of cool giants at resolving 
powers of at least 2000 (like the highest non-Fabry-Perot modes 
of ISOSWS and comparable to several groundbased instruments). 
Another is to rely partially on normal A-dwarfs whose spectra can 
be well modeled. Yet another is to extend the current network 
down from a 12-pm flux density of 5 Jy to 5 mJy. I illustrate these 
new directions through use of ISO spectra, KAO spectra of aster- 
oids, and the calibration needs of WIRE. 

References: [1] Cohen M. et al. (1999) Astrophys. J. t in 
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FIRST DETECTION OF POLARIZATION IN THE CO 
BANDHEAD EMISSION IN YOUNG STELLAR OBJECTS: 
A SIGNATURE OF THE INNER REGIONS OF DISKS IN 
SSV 13 AND DG TAU. C. Colome, Instituto de Astronomia, 
Universidad Nacional Autonoma de Mexico, Mexico, and 
McDonald Observatory, University of Texas- Austin, Austin TX, 
USA. 

We report the first detection of polarization in the CO band- 
head emission at near-infrared wavelengths (I = 2.3 mm). In this 
paper we describe the results for two young stellar objects: SSV 
13 and DG Tau. The near-infrared polarization (P > 13%) of the 
CO bandhead in both objects is considerably higher than the 
polarization in their continua. This suggests that, at least in those 
two objects, the CO emission originates in disks, rather than in 
stellar winds. 

It has been 20 years since the first detection of the first over- 
tone bands of CO in the Becklin-Neugebauer (BN) object [1]. Pos- 
sible mechanisms for the excitation of the vibrational states of CO 
have been described by Scoville et al. [2] and Krotkov et al. [3]. 
The detection of the CO bandhead in the BN object and further 
NIR spectroscopy of the same region [4] showed that in order to 
reproduce the relative band strengths, high vibrational temperatures 
are required (-3500 K). It was also noticed since then [1,4] that 
the CO bandhead emission must be localized within a small radius 
(a few AU) of the BN object, since the features are seen with the 
same strength in spectra taken with quite different apertures. 


TABLE 1. CO and continuum polarization. 


Object 

Filter 

P(±2%) 

Q(±3°) 

SSV13 

CO bandhead (v = 2-0) 

16 

66 


Kw’ (1. 9-2.4 mm) 

5 

73 

DG Tau 

CO bandhead (v = 2-0) 

13 

79 


Kw' (1.9-2 .4 mm) 

6 

108 


Dozens of young stellar objects (YSOs) have been detected via 
their CO bandhead, either in emission or absorption [5-11]. All 
the observations exhibit band strengths that are consistent with ex- 
citation temperatures of 2500-4000 K. It is still not known how 
this gas is excited and where it is located, nor it has been possible 
to study its kinematics. Two possible scenarios have been pro- 
posed for the origin of the CO emission: (1) an outflowing wind, 
or (2) a neutral circumstellar disk [7 and references therein]. Re- 
cent high-resolution infrared spectroscopy in the CO bandhead 
[10-12] of YSOs shows more evidence for a “disk” origin for the 
CO emission, rather than a wind origin. All these recent observa- 
tions show that the bandhead often exhibits the characteristic shape 
of emission from a rotating disk [12 and references therein]. A 
contribution from a wind to the CO emission of a rotating disk can 
be ruled out by the absence of line asymmetries in the NIR spec- 
tra. Line asymmetries would be expected due to the occultation the 
receding flow by an optically thick disk. An additional line asym- 
metry would be expected if the wind (cooler than the star and 
disk) would absorb the stellar and disk continuum. These effects 
have not been observed to date. 

The motivation for this work is the increasingly common de- 
tection of the CO bandhead emission in YSOs [7] since it was first 
discovered by Scoville et al. [1]. Most YSOs are surrounded by 
dust envelopes, which scatter not only the continua of young stars 
but also the radiation originating in the inner regions of their disks 
(if present). Scatter radiation can be detected as “polarized light.” 
Narrow-band polarimetry and centered at 2.3 mm, the first over- 
tone (v = 2-0) of the CO emission, of YSOs could provide per- 
haps the best method to distinguish between the two possible 
origins: circumstellar disks or neutral stellar (or disks) winds. 

As part of an ongoing observational project at McDonald 
Observatory SSV 13 DG Tau have been observed at near-infrared 
wavelengths using the polanmetnc near-infrared camera ROKCAM 
[13] of the University of Texas. The CO bandhead emission from 
both SSV 13 and DG Tau is unresolved with our resolution. Their 
FWHM are approximately 1.2" (3 pixels). For both objects the 
polarization measured with the narrow band filter centered in the 
CO bandhead is much higher than the one obtained using a 
broader K-band filter. The results are summarized in Table 1. For 
these two objects the polarimetry provides evidence for a “disk” 
origin for the CO emitting region, rather than a wind origin. In the 
latter case we would expect the scattering dust to see a more “iso- 
tropic” radiation from the CO gas, locally produced, and in such 
case, the large values of the polarization (P > 13%) are difficult 
to reproduce. 
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THE THERMAL INFRARED SPECTRA OF COMETS 
HALE-BOPP AND 103P/HARTLEY 2 OBSERVED WITH 
THE INFRARED SPACE OBSERVATORY. J Crovisier 1 , K 
Leech 2 , D. Bockelee-Morvan 1 , E. Lellouch 1 , T. Y. Brooke 3 , M. S. 
Hanner 3 , B. Altieri 2 , H. U. Keller 4 , T. Lim 2 < 5 , T. Encrenaz 1 , A. 
Salama 2 , M. Griffin 5 , T. de Graauw 6 , E. van Dishoeck 7 , and R. F. 
Knacke 8 , ^bservatoire de Paris-Meudon, F-92195 Meudon, 
France (crovisie@ obspm.fr), 2 Infrared Space Observatory Data 
Centre, Astrophysics Division of ESA, Villafranca, Spain, 3 Jet 
Propulsion Laboratory, Pasadena CA, USA, 4 Max-Planck-Institut 
fur Aeronomie, Katlenburg-Lindau, Germany, 5 Queen Mary and 
Westfield College, London, UK, 6 SRON, Groningen, The 
Netherlands, 7 Leiden Observatory, The Netherlands, 8 Pennsylvama 
State University, Erie PN, USA. 

Introduction: The Infrared Space Observatory (ISO) offered 
us the opportunity to observe celestial bodies over the full infra- 
red range from 2.4 to 196 pm. This spectral domain is of pecu- 
liar interest for comets. Cometary parent molecules, likely to have 
sublimated from the nucleus ices, can be investigated through their 
fundamental bands of vibration in the 2.5-12-pm region. With 
temperatures ranging between 100 and 400 K or more, cometary 
grains show thermal emission peaking in the 5-200-pm region, 
where spectral features related to their chemical nature may be 
expected. 

Observations: We report here on spectroscopic observations 

of comets that were conducted with ISO by an international team 
of investigators. Other programs by other teams were aimed at 
photometric and imaging observations, in order to study cometary 
dust in the coma, the tails and the trails, and to try to disentangle 
the emission of the comet nuclei from that of the dust. 

Our observations were undertaken in the frame of two pro- 
grams. The first one, a guaranteed time program, had to be 
planned on advance. It was thus aimed at short-period, predictable 
comets. The nominal target, 22P/Kopff. was weaker than expected 
and could not be observed at the moment it was the brightest. 
Fortunately, it was possible to observe a second comet, 103P/ 
Hartley 2, during the extended lifetime of ISO, when the comet 
was at its brightest near perihelion [7]. The second program, an 
open-time target-of-opportunity program, was aimed at any bright 
comet that might appear during ISO’s lifetime. We were lucky to 
benefit from the apparition of the exceptional Comet Hale-Bopp 
(1995 01). Unfortunately, because of the visibility constraints of 
ISO, the comet could only be observed when it was at r h > 2.8 AU 
from the Sun. Observations were secured preperihelion on April 
1996 (r h = 4.6 AU) and September -October 1996 (r h = 2.8 AU) 
and postperihelion on December 1997 (r h = 3.9 AU) and April 
1998 (r h = 4.9 AU) [1-6]. 

The observations of gas-phase species will not be discussed in 
detail here: H-,0, C0 2 , and CO were detected, and their produc- 


tion rates, as well as physical parameters for cometary water, were 
determined as a function of heliocentric distance [1— 4,6,7]. 

Results on Comet Hale-Bopp: For the four observing peri- 

ods of Comet Hale-Bopp, low-resolution spectra were observed 
from 2.5 to 12 jam with the grating spectrometer of the ISO pho- 
tometer. All spectra show a strong silicate band at 10 pm with a 
peak at 11.3 pm indicating the presence of crystalline olivine. This 
constituent is observed to be present whatever the heliocentric 
distance (4.9 to 2.8 AU), both pre- and postpenhelion. 

In September-October 1996, high signal -to-noise ratio spectra 
of the comet were obtained with the short- and long-wavelength 
spectrometers of ISO. Beyond 6 pm, the spectrum of Comet Hale- 
Bopp was dominated by dust thermal continuum emission, upon 
which broad emission features were superimposed, the strongest 
appearing at 10, 19.5, 23.5, and 33.5 pm. The wavelengths of all 
these peaks correspond to those of Mg-rich crystalline olivine 
(namely forsterite, Mg 2 Si0 4 ) [2-4,6]. A careful study of the ISO 
spectrum, complemented by groundbased spectra of the 10-pm 
region obtained at smaller heliocentric distances, show that amor- 
phous silicates as well as pyroxenes are also present [8,9]. It is 
remarkable that the spectrum of Comet Hale-Bopp in this spectral 
domain closely resembles that of the dust disks around some 
young stars (such as that of HD 100546, an intermediate star bet- 
ween Vega-like objects and Herbig Ae-Be stars) that were also ob- 
served by ISO [10]. Thus, the dust contained in such disks should 
be similar to cometary dust. This contrasts with interstellar sili- 
cates, which are amorphous. 

In the September-October 1996 spectra, emission features at 
44 and 65 pm and possible absorption at 2. 9-3. 2 pm were detec- 
ted, characteristic of crystalline water ice. A temperature of 165 K 
was inferred for this water ice. This suggests that grains of water 
ice were still present at a distance of 2.9 AU from the Sun and that 
they could be a significant source of water release in the coma [5]. 

No sign of PAHs could be found in the spectrum of Comet 
Hale-Bopp (especially in the 6-9 pm region, where they are ex- 
pected to have strong bands not blended with silicate features). 
However, one cannot exclude that cometary PAHs are revealed at 
shorter heliocentric distances than those of the ISO observations 
of Comet Hale-Bopp. 

Results on Comet P/Hartley 2: Whereas Hale-Bopp was a 

long-period comet presumably coming from the Oort cloud, 103P/ 
Hartley 2 is a Jupiter-family comet that supposedly originated from 
the Edgeworth-Kuiper belt. It was observed in December 1997— 
January 1998 at r h = 1 AU. A 6-17-pm low-resolution spectrum 
obtained with the VF of ISOCAM showed a weak 10-pm silicate 
band with an intensity only 15% of the continuum [7]. The pre- 
sence of a peak at 11.3 pm revealed the presence of crystalline 
silicates 

Conclusion: ISO observations — complemented by ground- 

based spectroscopic observations in the infrared and radio domains 
— show that there is a great similarity in the composition of in- 
terstellar and cometary ices. This agrees with the hypothesis that 
comets could have formed from basically unprocessed interstellar 
matter. However, the presence of crystalline silicates in the dust 
of both a Oort-cloud comet and a Jupiter- family comet shows that 
these bodies have also incorporated processed material at some 
stage of their formation. 
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AIRES: AN AIRBORNE INFRARED ECHELLE SPECTRO- 
METER FOR SOFIA. J. Dotson 1 , E. Erickson 2 , M. Haas 2 , S. 
Colgan 2 , J. Simpson 1 , C. Telesco 3 , J. Wolf 4 , and E. Young 5 , 
! SETI Institute, Mail Stop 245-6, NASA Ames Research Center, 
Moffett Field CA 94035, USA (dotson@cygnus.arc.nasa.gov), 
2 NASA Ames Research Center, Moffett Field CA 94035, USA, 
3 Department of Astronomy, Bryant Space Science Building, Uni- 
versity of Florida, Gainsville FL 32611, USA, 4 DLR, Germany, 
5 University of Arizona, Tucson AZ, USA. 

SOFIA will enable astronomical observations with unprec- 
edented angular resolution at infrared wavelengths obscured from 
the ground. To help open this new’ chapter in the exploration of 
the infrared universe, we are building AIRES, an Airborne Infra- 
Red Echelle Spectrometer. AIRES will be operated as a first gen- 
eration, general purpose facility instrument by USRA, NASA’s 
prime contractor for SOFIA. AIRES is a long-slit spectrograph 
operating from 17-210 pm. In high-resolution mode the spectral 
resolving power is -10 6 pmA or -10 4 at 100 pm. In low-resolu- 
tion mode the resolving power is about lOOx lower. AIRES in- 
cludes a slit-viewing camera that operates in broad bands at 18 and 
25 pm. 

AERES will be ideal for spectral imaging of gas-phase phenom- 
ena in the interstellar medium (ISM). Far-infrared line observations 
probe the excitation, pressure, density, luminosity, chemical com- 
position, heating and cooling rates, mass distribution, and kine- 
matics in the various components of the ISM. The lines offer 
invaluable and often unique diagnostics of conditions in such 
diverse places as star-forming regions, circumstellar shells, the 
Galactic Center, starbursts in galaxies, and the nuclei of active 
galaxies. AIRES will provide astronomers with new insights into 
these and other environments in the ISM. It will also be useful for 
studies of solar system phenomena such as planetary atmospheres 
and comets, for example. 

AIRES’ low-resolution mode will enable comparison of infra- 
red continuum emission with line emission from interstellar clouds. 
It will also permit measurement of long wavelength features of 
interstellar grains and solid surfaces of solar system bodies, as well 
as a variety of other problems. 
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AXISYMMETRIC MODELING OF SILICATE CIRCUM- 
STELLAR DISKS. S. K. Dunkin and M. J. Barlow', Depart- 
ment of Physics and Astronomy, University College London, 
Gower Street, London WC1E 6BT, UK (skd@star.ucl.ac.uk). 

Introduction: The presence of dust grains around many stars 

is not disputed. However, the composition and size of these dust 
grains, as well as the geometry and inclination the dust takes as 


a whole, is. There are three main classes of infrared (IR) excess 
star whose dust may be representative of either the beginning or 
remnants of planetary systems: the main- sequence Vega-like stars 
and pre-main sequence Herbig Ae/Be (HAeBe) and T Tauri stars. 

There have been many attempts to model the spectral energy 
distributions (SEDs) of HAeBe stars and their lower mass coun- 
terparts to determine the characteristics of the dust [e.g., 1-4]. 
Attention has also been paid to the prototype Vega- like stars too, 
which exhibit only a mid- or far-IR excess, but relatively little 
work has been carried out on the class of Vega-like stars that have 
substantial near-IR excess as well. The excess emission from some 
of these stars cannot be modeled completely using optically thin 
codes [5,6], and this work concentrates on modeling one of them 
using an optically thick radiative transfer code. 

Modeling: Several Vega-like stars from the list of [7] were 

studied and one, HD 144432, with both a near- and far-IR excess, 
was chosen for modeling. This star has a very strong silicate 
emission feature (Fig. 1 and [8]). Such a signature helps to place 
strong constraints on the compositional parameter of the model. 

The code used here is a modified version of that used by [9]. 
Skinner et al. [10] explains the major refinements of the code, and 
Dunkin [11] describes in more detail the various parameters and 
usage of it. There are many free parameters in this code, but it is 
still possible to place constraints on certain values. 

Spherical vs. nonspherical distributions . While asymmetric 
geometries have been positively identified for a few stars (e.g., p 
Pic [12], HD 163296, and AB Aur [13]), the dust geometries of 
most IR excess stars remains uncertain. Before attempting model- 
ling, it is possible to provide evidence for or against the presence 
of spherical envelopes by comparing the extinction observed to- 
wards a star to that expected if such an envelope was present. By 
using this method (as outlined in [13]) we determined that only 
two of the 1 1 Vega-like stars we looked at could potentially have 
spherical distributions of dust, neither of which was our modeled 
star, HD 144432. Miroshnichenko et al. [4] have been able to pro- 
duce fits to several HAeBe stars using spherical envelopes with 



Fig. 1. Optically thick silicate model of HD 144432 plotted on observa- 
tional data. Note the large silicate feature at 10 pm. 
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dust density profiles obeying an r* 1 5 power law. Their model is 
sufficiently similar to the one used here for comparisons to be 
made between the two. We created models for two of their stars 
using parameters as close as possible to theirs. We confirmed that 
spherical envelopes may produce reasonable fits to the data for 
these stars, and used these parameters as a base for the modeling 
of our chosen star. No spherical model could be made to fit the 
observed data for HD 144432, despite changing inclination, grain 
size distribution, optical depth, density distribution, and composi- 
tion. We therefore conclude that HD 144432 has its dust distrib- 
uted in a nonspherical geometry, most likely in a disk. 

Composition of the disk. Sylvester et al. showed HD 144432 
to have a very strong silicate feature [8]. Models were run using 
composite grains of amorphous C and silicate, but only a model 
using 100% silicate grains could reproduce the observed silicate 
peak. It is unlikely that the actual disk is pure silicate, but this 
model does show that silicate grains must dominate the total 
population. 

Other parameters . A grain size distribution was used as op- 
posed to a single grain size as this is probably a more realistic 
situation. Interestingly, millimeter-sized grains were found to be 
necessary to produce the flux seen at millimeter wavelengths. This 
implies that at least some of the dust has accreted to form large 
grains. We cannot put an upper limit on the grain growth in this 
system, but cannot rule out the possibility of even larger grains 
being present. An inner radius of the disk at 0.13 AU was required 
to produce the near-IR excess: If the inner regions of the disk 
were further out, not enough flux at these short wavelengths can 
be produced. Farther out, large grains and an outer radius of 
-1300 AU is required to reproduce the millimeter-wave flux. The 
outer radius is more uncertain than the inner radius, and may well 
be larger. 1300 AU in our solar system places us somewhere in 
or just outside the Kuiper belt where a lot of “left-over” debris 
from the formation of our planets resides. However, HD 144432 
is an A9/F0 star and is therefore hotter and more massive than our 
own Sun, so direct comparisons between systems is difficult to 
make. The disk is optically thick, with an optical depth, x, of 3.0 
required to produce the flux seen. The disk is therefore quite 
substantial, with a mass of 2.71 x 10" 6 solar masses. Since the 
inclination of the system is unknown, an inclination of 45° was 
assumed (edge-on and face-on inclinations are ruled out by the 
vsim of the star). Inclination is known to have little effect on the 
resultant models using this code, unless it is close to 0° or 90° 
[11]. The final model is shown in Fig. 1 below. 

Summary: In summary, we have used an optically thick ra- 

diative transfer code to model a Vega-like star, HD 144432, with 
near-IR excess. The dust around this star is not distributed in a 
spherical geometry and can be fit well using a flattened disk dis- 
tribution. The large silicate feature at 10 pm implies that the dust 
in this system is dominated by silicate grains, a fact reflected in 
the final model. 

This code can best be used in situations where the dust is likely 
to be optically thick, as in the case of HD 144432. Further work 
is underway to model other stars, both main sequence and pre- 
main sequence, to compare with the results presented here. 
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LABORATORY MEASUREMENTS OF SPECTRAL EMIS- 
SIVITY OF MATERIALS ANALOGS TO THE MARTIAN 
SOIL AND DUST, S. Fonti 1 , A. Blanco 1 , M. I. Blecka 2 , F. 
De Carlo 1 , V. Orofino 1 , and S. Solazzo 1 , Physics Department, 
University of Lecce, C.P. 193, 73100 Lecce, Italy (fonti@le. 
infn.it), 2 Space Research Centre, Polish Academy of Sciences, ul. 
Bartycka 18a, 00-176 Warsaw, Poland. 

Introduction: In recent years the debate about the climatic 

history of Mars has been quite alive. In particular the presence of 
liquid water in large basins for long periods is supported by sev- 
eral indications [1,2]. One of the main issues is to understand if 
such climatic conditions lasted long enough to allow the develop- 
ment of some form of life [3]. Although it seems premature to 
expect a definite answer on the subject, interesting suggestions 
could be given by a systematic spectroscopic survey of the mar- 
tian surface. As an example, the detection of carbonates in selected 
regions of the planet [4], can be important to assess that the pres- 
ence of large bodies of standing water has lasted long enough to 
allow an abundant precipitation of such compounds. 

Unfortunately, up to now no definite identification has been 
made by any instrument orbiting around Mars. In fact, the Ther- 
mal Emission Spectrometer (TES), onboard Mars Global Surveyor, 
has set only an upper limit of about 10% to the amount of these 
materials [5], while the detection, yet to be confirmed, about the 
possible presence of carbonates comes from two observations 
made from the Earth [6,7]. In these cases, however, it is almost 
impossible to link the observation to a specific area of the planet. 

In this scenario an important role can be played by laboratory 
measurements on materials, such as carbonates, which are sup- 
posed to be present on the martian surface, as well as in the at- 
mospheric dust. The measurements should be carried out in 
different environmental conditions, on both bulk and dust samples, 
varying not only the composition, but also the average grain size 
and the compactness of the granulated specimen [8]. The output 
of the laboratory measurements can be fed into an appropriate 
radiation transfer model able to produce reliable synthetic spectra 
of the martian surface and atmosphere. Such spectra could be 
conveniently compared with the actual spectra of Mars, which will 
be available from current (Mars Global Surveyor) and future (Mars 
Express) martian missions for the identification of the spectro- 
scopic signatures of materials. In this respect it is worthwhile to 
note that most of such materials (carbonates, sulfates, silicates, 
etc.) exhibit their characteristic bands in the mid-infrared region 
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of the spectrum (5-50 pm) and that in this range the martian 
spectra are dominated by the thermal emission of the surface. 

The emissivity measurements presented in this work are part of 
a research program, carried out in the Astrophysics Laboratory of 
the Department of Physics of the University of Lecce, aimed to a 
systematic study of the possible analogs to the martian materials. 

Experiment: The measurements of spectral emissivity have 
been performed using a Perkin-Elmer Fourier Spectrometer Model 
Spectrum 2000, equipped with an emission port. The spectral 
range was limited between 2000 and 500 cm* 1 (5-20 pm), due to 
the presence of high noise level outside this interval. It would be 
useful to extend the range on the long wavelength side and we are 
presently working in this direction. An optomechanical accessory 
was especially designed and manufactured in our laboratory, in 
order to concentrate the energy emitted by the sample into the 
input port of the spectrometer. Such accessory allows, by means 
of a moveable mirror, to look alternatively to the sample, kept at 
a nominal temperature variable between ambient and 300°C, a 
reference black body at the same nominal temperature, and a sec- 
ond black body at constant ambient temperature. 

The comparison between the spectrum obtained looking at the 
sample and that of the hot black body is the basic process for 
calculating the spectral emissivity, while the spectrum of the cold 
black body allows to account for the real zero of the spectrometer. 
However, the actual procedure is more complicated and requires 
several intermediate steps in order to assess some important pre- 
liminary parameters. As an example, we can consider the unavoid- 
able difference between the nominal temperature at which the 
thermostatic devices are set and the real temperature of the emit- 
ting surfaces. To account for this, it is necessary to perform sev- 
eral measurements at different nominal temperature and compare 
the results. The effective temperature of the emitting surfaces is 
found by means of an iterative calculation, based on the compari- 
son between the experimental results and the theoretical predic- 
tions. 
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Fig. 1. Synthetic spectra of Mars, calculated in vertical geometry, 
obtained using the laboratory measured emissivity of palagonite and 
limestone. The two spectra refer to a surface of pure limestone (thick line) 
and of pure palagonite (thin line); in both cases the atmosphere is thought 
to be loaded with micrometer- si zed grains of pure palagonite. 


It is important to stress that, although the temperature can be 
varied in a wide range, the best results have been obtained for a 
nominal temperature higher than 1 50°C. Even if other authors [9- 
1 1] have done emissivity measurements using a considerably lower 
temperature, our choice has the great advantage to reduce drasti- 
cally the influence of the reflected radiation on the result of the 
measurement. On the other hand, using a temperature much higher 
than the ambient temperature can produce temporal instability and 
spatial inhomogeneity. For the second point, the use of the effec- 
tive temperature of the emitting surface, directly evaluated from 
the measurements, will conveniently average on possible small 
inhomogeneities. As far as the temporal instability is concerned, 
we think that the use of a fast instrument, such as a Fourier spec- 
trometer, together with a continuous monitoring of the black body 
kept at ambient temperature, ensures that the level of a possible 
temporal instability will be within a few percent. 

Results and Discussion: In order to perform a preliminary 

check of the visibility of the carbonate component in the spectra 
of Mars, we measured the spectral emissivity of a typical carbon- 
ate as limestone and of palagonite, which is assumed to be the best 
terrestrial analog to martian silicates. We then fed our measured 
emissivities in an advanced model of radiative transfer calculation 
[12], which takes into account the emission of the martian surface 
and the emission and the extinction by the atmospheric dust and 
gases. The resulting synthetic spectra, shown in Fig. 1, have been 
obtained using as main input parameters to the model the vertical 
optical thickness of the dust at 10 pm (x D = 1.1), the pressure 
(p s = 6 mbar), and the temperature (T s = 270 K) at the surface of 
Mars. A more detailed discussion of the model can be found else- 
where [13). 

As it can be seen, two main features characterize the spectrum 
obtained assuming a pure carbonate surface (thick line): the typi- 
cal band at 6.7 pm and another minor, but still evident structure 
at 1 1 .4 pm. Of course, these features are not seen in the calcula- 
tion assuming pure palagonite surface (thin line). It is interesting 
to note that these features are detectable even when the atmosphere 
is loaded with palagonite dust. Obviously, in order to draw a 
specific conclusion, more tests have to be done, using composite 
surfaces with variable percentages of the main components. We 
can conclude, however, that the results of the present work indi- 
cate that the comparison of the synthetic spectra with the obser- 
vational ones can be a very promising tool for the spectroscopic 
identification of carbonates and other materials on the surface of 
Mars. 

Acknowledgments: This work has been partially supported 

by the Italian Space Agency and by the Polish State Committee for 
Scientific Research under grants 2P03C01010 and P03C00515. 

References: [1] Goldspiel J. M. and Squyres S. W. (1991) 

Icarus, 89, 392-410. [2] Scott D. H. et al. (1995) USGS Misc. 
Inv. Ser . Map 1-246 L [3] McKay D. S. et al. (1996) Science, 273, 
924-930. [4] Orofino V. et al. (1999) Proceedings of the 2nd 
Italian Meeting of Planetary Sciences, Bormio, 20-25 January 
1998, in press. [5] Christensen P. R. et al. (1998) Science, 279, 
1692-1698. [6] Pollack J. B. et al. (1990) JGR, 95, 14595-14627. 
[7] Leilouch E. et al. (1999) Proceedings of the International 
Symposium on Mars Exploration Program and Sample Return 
Mission, Paris, 2-5 February 1999, 011/52 (40). [8] Arnold G. 
(1991) Vibrational Spectroscopy, 2, 245-249. [9] Nerry F. et al. 
(1990) JGR, 95, 7027-7044. [10] Christensen P. R. and Thlivens 




LPl Contribution No. 969 11 


Harrison S. (1993) JGR, 98, 19819-19834. [11] Salisbury J. W. 
et al. (1994) JGR, 99, 11897-11911. [12] Blecka M 1. et al. 
(1997) Adv. Space Res., 19, 1281-1284. [13] Blecka M. 1. (1999) 
Adv. Space Res., in press. 


THEORETICAL MODELING OF RADIATIVE TRANSFER 
IN PLANETARY REGOLITHS. B. W. Hapke, Department of 
Geology and Planetary Science, University of Pittsburgh, Pitts- 
burgh PA 15260, USA. 

The status of theories of thermal emission from particulate 
media, such as planetary regoliths, is discussed. When the particles 
are large compared with the wavelength, models based on the 
radiative transfer equation (RTE), which assumes that the particles 
scatter independently, appear adequate. Geometric optics can be 
used to calculate the optical properties of the particles because 
diffraction by the particles effectively disappears. The emissivities 
of individual particles can display spectral features of high con- 
trast. However, multiple scattering in a regolith may greatly reduce 
the contrast, decreasing the usefulness of thermal IR spectral 
measurements for remote compositional identifications [1]. The 
maximum in the spectral emissivity does not occur exactly at the 
Christiansen frequency, as is commonly assumed (although it is 
close), but rather at the wavelength in a reststrahlen band where 
the light scattered from a particle changes over from scattering 
within the interior of the particle (volume scattering) to reflection 
from the surface [2]. 

In order to calculate the emissivity and subsurface temperature, 
the RTE must be simultaneously solved along with the heat trans- 
fer equation [3]. Such solutions intrinsically include radiative 
conductivity and the solid state greenhouse effect. Including radia- 
tive effects quantitatively alters the calculated subsurface tempera- 
ture distribution. 

The status of thermal emission models is much less satisfactory 
when the particle size is comparable to or smaller than the wave- 
length. The use of the RTE is questionable then because coherent 
effects do not allow the particles to scatter independently. Unfor- 
tunately, this case probably includes most planetary regoliths. 
Since coherent effects have a range of roughly a wavelength, this 
suggests an approach in which the regolith is averaged over wave- 
length-sized boxes to obtain effective particles whose optical 
peoperties may be inserted into the RTE. 
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OPTICAL PROPERTIES OF COSMIC DUST ANALOGS. 

Th. Henning, Astrophysical Institute and University Observatory, 
Friedrich Schiller University, Schillergasschen 3, D-07745 Jena, 
Germany (henning@astro.uni-jena.de). 

Small solid particles form an important component of our 
galaxy and others. Although they are only minor contributors to 
the mass of these systems, they determine the thermal, dynamical, 
and chemical state of the dense phases of the interstellar medium, 
thereby controlling the star formation process. The coagulation of 


these particles in circumstellar disks is a first but decisive step 
toward the formation of planetesimals. 

Cosmic dust grains are a very interesting system of small par- 
ticles interacting with radiation, other grains, and gas particles 
under the extreme conditions of interstellar space. Laboratory 
investigations have to supply the basic data for these interactions 
in order to interprete the results of ultraviolet, optical, and infra- 
red spectroscopy of dusty objects. 

The talk will review laboratory measurements of the optical 
data of cosmic dust analogs. It will give an introduction how such 
materials are produced and characterized. The main techniques to 
measure optical constants with their pros and cons will be dis- 
cussed. 

Computational methods to calculate absorption and scattering 
efficiencies for fluffy and inhomogeneous particles from the op- 
tical constants will be presented. 
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PROBLEMS IN INFRARED LABORATORY SPECTROS- 
COPY OF SILICON CARBIDE RESULTING FROM HIGH 
EXTINCTION COEFFICIENTS: INFERENCES FROM 
NANOCRYSTALS AND QUANTITATIVE METHODS. 

A. M. Hofmeister, 1 A. K. Speck 2 , and L. Rosen 3 , Apartment of 
Earth and Planetary Science, Washington University, St. Louis 
MO, USA (hofmeist@levee.wustl.edu), ^Department of Physics 
and Astronomy, University College London, UK, 3 Department of 
Mechanical Engineering, Washington University, St. Louis MO, 
USA. 

Introduction: Infrared measurements of minerals in the labo- 

ratory are needed for identification of phases detected remotely in 
comets or in circumstellar dust. Silicon carbide is of particular in- 
terest due to its association with cool, C-rich stars. Unfortunately, 
problems with convolution of the baseline spectrum with that of 
the sample have gone unrecognized. As a result, spectral artifacts 
have been attributed to intrinsic features. Such problems are ac- 
centuated for SiC due to high reflectivity [1] and consequent high 
absorptivity [2]. 

It is generally assumed that micrometer or submicrometer par- 
ticles transmit intense IR bands. As a consequence, the weight of 
the particles in the pellet have been used to determine the extinc- 
tion coefficient, and the extinction coefficients inferred from 
single-crystal reflectance studies have been discounted because 
these were too large in comparison. However, recent study of 
single crystals of grossular and andradite garnets thinned by ion- 
milling and of powders compressed into a thin film in a diamond 
anvil cell showed that thicknesses below 0.1 pm are needed to 
transmit the most intense Si-0 modes. Moreover, the results are 
consistent with absorptivity calculated from reflectance spectra 
[3,4]. These data indicate that intense peaks in a KBr pellet are 
due only to the finest grains, or arise from edges of grains. Be- 
cause “opaque” grains exist in the sample, the baseline depends on 
the sample spectrum, and as a result, peak positions, widths, and 
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shape are falsified when the ratio of the raw data to a reference 
pellet is calculated. 

Related problems exist in IR spectra from small particles ob- 
tained with an IR microscope. Because the image is typically 
masked, rather than the sample, light leakage (e.g., diffraction) 
around the edges of essentially opaque particles alters the zero 
percent transmission line [5,6]. The resulting spectra are problem- 
atic. 

Experiments: To further our understanding of extinction 

coefficients and to delineate circumstances under which spectra are 
intrinsic, rather than circumstantial, we examined IR spectra from 
small single crystals of a-SiC, nanocrystals of structurally related 
AIN and TiN, and (3-SiC. and commercial polycrystals. All nano- 
samples have grains with 2-10-nm sizes. 

The nanocrystals of AIN and TiN were encapsulated in NaCl 
during crystallization. Standard pellets produced very broad bands, 
however, if the sample was diluted by compression, the spectra 
gradually changed into a moderately sharp peak, the position of 
which was replicated by a thin film of commercial, bulk sample. 
The production of same spectra of nano and bulk sample confirms 
the high extinction coefficients, the production of IR modes essen- 
tially by nearest neighbor interactions, the inconsequential contri- 
bution of surface modes, and the neutral role of the salt matrix. 

Nanocrystals of p-SiC occur as clusters without a salt matrix. 
Grinding reduced the cluster size to micrometers as measured by 
TEM of one of our KBr pellets. Their thin film spectra are domi- 
nated by the TO mode at 12.2 pm, and dispersion spectra have TO 
and LO modes subequal in intensity. The dispersion spectrum 
most closely resembles thin film spectra of commercial p sample 
[2]. These results are consistent with the extreme hardness of SiC 
and initial crystal sizes. For the nanocrystals, grinding did not 
reduce the cluster size to sufficiently small sizes to obtain an 
intrinsic spectrum, but crushing did, whereas macroscopic crystals 
of p-SiC could not be sufficiently thinned even under pressure. 
Further experiments with smaller clusters are underway. 

Hexagonal a-SiC platelets, 150 jim across by 5 pm thick, have 
maximum absorbance >2 units when sample is placed on a 
100-pm aperture, but the same cystal yields a flat topped peak at 
0.5 units if the image is masked to an equivalent aperture instead. 
For 200-pm-diameter cystals the disparity is less, but still exists. 
This suggests that diameters of the sample must be >4x that of the 
effective aperture for spectra to be representative of the sample. 

Reflectivity measurements of SiC repeated for E perpendicular 
to c axis confirm the results of Spitzer et al. [1]. Their absorbance 
data does not agree with our above measurements, probably be- 
cause their film was spotty. However, their calculate extinction 
coefficients are confirmed. 

Discussion: All the above data, taken together, indicate that 

for SiC particles to transmit all wavelengths, the path length must 
be considerably below 0.1 pm. For larger path lengths, peak height 
and profiles are altered through instrumental errors, specifically by 
convolution of the baseline with the spectrum. For thickness above 
a micrometer, the spectra are essentially all artifacts. Hence, most 
published absorbance spectra of SiC have artifacts to varying 
degrees, and connection of trends in peak parameters with factors 
such as grain size, grain shape, and use of a matrix are unsup- 
ported. 

Reasonable data can be obtained from dispersions if grain size 
distribution is restricted to small sizes and a narrow size range, as 


well as to dilute samples, as in the finest particle distributions 
published to date [7]. Use of heavy dispersions and maximum 
peak heights above 0.5 to 1 absorbance units distort the peak 
profiles because the total path length is likely to be opaque at the 
peak center, while still transparant to the shoulders or sides of the 
peak. 

It is recommended to use a thin film [2], as path length is 
controlled, and the same sample can be thinned until the spectrum 
no longer changes. Typically, the spectra no longer change once 
peak height is below 0.5 absorbance units. This effect occurs 
because the thin films are not perfectly uniform. 

Conclusions: Comparing laboratory spectra to circumstellar 
spectra requires caution as the baseline subtraction procedures 
differ. However, if both the absorptions from the laboratory and 
in astonomical data are weak, the problems should be minimized. 
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TRANSIENT DUST FEATURES IN CARBON-RICH CIR- 
CUMSTELLAR ENVELOPES OF EVOLVED STARS. S 

Hony 1 , R. J. Sylvester 2 , L. B. F. M. Waters 1 * 3 , and M. J. Barlow 2 , 
Astronomical Institute, University of Amsterdam, Kruislaan 403, 
1098 SJ Amsterdam, The Netherlands (hony@astro.uva.nl), 
department of Physics and Astronomy, University College 
London, London WC1E 6BT, UK, 3 SRON Laboratory for Space 
Research Groningen, P.O. Box 800, 9700 AV Groningen, The 
Netherlands. 

Introduction: The last stages in the life cycle of intermedi- 

ate mass stars, with masses ranging between 1 and lOx the mass 
of our Sun is characterized by a very rapid evolution in structure, 
luminosity, and surface temperature. Changes in luminosity of the 
asymptotic giant branch stars (AGB stars) induce convection in the 
outer parts of the stars, causing nuclear burning products from the 
CNO cycle to be transported to the surface. This causes the abun- 
dances of the various elements in the atmosphere of these stars to 
change. When, by this process, the ratio of C- over O-atoms be- 
comes larger then unity it changes the chemistry dramatically and 
the atmosphere becomes C-rich. This is caused by the stability of 
the CO molecule, which locks up most of the C- and O-atoms, and 
only the more abundant species is available to form other mol- 
ecules. 

At the same time these stars are affected by strong mass loss, 
sometimes reaching values of 10 _4 @/yr. In the resulting high-den- 
sity outflows dust will condensate into grains with a composition 
dependent on the abundances of the available species. 

From AGB Stars to Planetary Nebula: When the outflow 

stops the gas and dust move away from the star, the star contracts 
and heats up from T - 10 3 K to T - 10 5 K with a timescale of 
-10,000 yr. When the star becomes hot enough, T - 10 4 K, to 
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ionize the expelled gas and dust, the system appears as a planetary 
nebula. During this rapid evolution several transient solid-state 
emission bands are observed in the infrared spectra from these 
objects. The main features observed are 6-9 and 10-15 pm attrib- 
uted to large PAH or HAC [1]; 21 pm unidentified [4]; and 24- 
50 pm unidentified, possibly MgS [3]. 

Transient Features: In this poster we study the properties of 

the dust in the rapidly changing physical conditions of these cir- 
cumstellar environments. We present new infrared observations, 
from AGB stars to planetary nebulae, secured with both the Short 
Wavelength Spectrometer (SWS) and the Long Wavelength Spec- 
trometer (LWS) onboard ESA’s Infrared Space Observatory (ISO). 
Here we focus on a broad solid-state emission feature ranging 
from 24 to 50 pm peaking around 30 pm. This 30-pm feature was 
first reported to be present in the spectra of C-rich stars by Forrest 
et al. [2]. Goebel and Moseley [3] suggested that solid MgS could 
be the carrier of this feature. However, because of the broadness 
of the feature and the difficulties of observing at these wave- 
lengths, no positive identification has yet been possible. 

The very large wavelength coverage of the two ISO instruments 
combined (2-200 pm) allows us to derive reliable band strengths 
and shapes for the 30-pm feature. We detect a change in shape 
and peak position of the feature with evolutionary status of the 
central star, leading us to conclude that at least two separate com- 
ponents are visible. The first feature is strongest in the AGB stars 
and peaks near 30 pm. The second component peaks near 35 pm 
and becomes more prominent in the planetary nebulae. We com- 
pare the observed spectra with model spectra of MgS and graph- 
ite grains and discuss the possibility of grain destruction by 
ultraviolet photons of the stellar radiation field and the effects of 
grain temperature on the emitted radiation. For identification pur- 
poses, we discuss the correlation of the these components with the 
other main solid-state emission features at shorter wavelengths. 
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SILICATES AS PROBES OF THE MASS LOSS HISTORY 
OF OXYGEN-RICH EVOLVED STARS. C Kemper 1 , R J 
Sylvester 2 , L. B. F. M. Waters 1 ’ 3 , T. de Jong 14 , M. J. Barlow- 2 , 
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Netherlands. 

Introduction: In the late stages of stellar evolution, interme- 

diate mass stars undergo severe mass loss (up to lO' 4 0 /yr) while 
ascending the asymptotic giant branch (AGB). The outer layers of 
the star are enriched with elements produced by the nuclear pro- 
cesses in the core of the evolved star. These elements include C, 
N, and O, where in most cases O is more abundant than C. When 
the circumstellar material cools due to expansion, gaseous SiO, 


together with Fe and Mg, will condense, thus forming silicates. As 
observations have shown [1], we can identify amorphous and 
crystalline silicates, whereas the different types probably originate 
at different temperatures and densities. 

One reason to study the infrared spectrum of AGB stars is that 
the emission and absorption features of crystalline and amorphous 
silicates will provide us with information on the physical and 
chemical conditions in the circumstellar dust shell. However, this 
is still future work. We intend to derive the temperature and den- 
sity gradients, thus revealing the mass loss history of O-rich AGB 
stars, and finally we may be able to put constraints on the nature 
of the mass loss mechanism. Here we will present spectra of O- 
rich evolved stars obtained with the Infrared Space Observatory 
(ISO). We will analyze the relation between the spectral features 
of silicates and the energy distribution. 

Energy Distribution: The dust shell surrounding the central 

star causes a redistribution of the emission. Radiation from the star 
in the visible wavelength region is (partly) absorbed by the dust 
and reemitted in the infrared region, thus causing an infrared 
excess. In case of a high mass loss rate and high column densi- 
ties toward the central star, the dust shell will be optically thick 
at visible wavelengths and sometimes even in the near-infrared. 
The central star is then completely obscured. This is the case for 
OH/IR stars. Bedijn [2] models the energy distribution using the 
dust shell surrounding O-rich stars due to (accelerated) mass loss. 
In this study the optical properties of amorphous silicate are taken 
into account. This provides a measure of the optical depth, since 
the amorphous silicate features at 9.7 and 18 pm are in emission 
for low column densities. The features go into self-absorption and 
absorption for increasing column density. Moreover, for increas- 
ing optical depth, the x = 1 surface will migrate outward to lower 
temperatures, which causes a redward shift of the peak of the 
energy distribution for increasing mass loss. 

From observations performed with ISO it becomes clear that 
the crystalline silicate features will become stronger for higher 
integrated mass loss. It is generally believed that annealing of the 
amorphous material into crystalline dust will happen only at high 
temperatures, i.e., above the glass temperature. In high mass loss 
objects, the density at the inner edge of the dust shell is higher. 
In that case, the dust condensation may start at a higher tempera- 
tures and the dust stays at those high temperatures for a longer 
period. Thus we increase the duration of the annealing period, 
which leads to a higher abundance of crystalline silicates. 

Observations: We present the combined Long Wavelength 

Spectrometer (LWS) and Short Wavelength Spectrometer (SWS) 
ISO data of Miras and OH/IR stars with different mass loss rates, 
e.g., OH26.5, WX Psc, OH104.9 and OH127.8. The infrared spec- 
tra span the wavelength range from 2.4 to 195 pm, which includes 
practically all the flux emitted by the sources. 

Analysis: The objects are ordered with increasing optical 

depth, using the 10-pm silicate band, as well as the peak position 
of the energy distribution. We find that the crystalline silicate 
features become stronger for cooler shells. The pseudocontinuum, 
i.e., the thermal emission from the dust shell, is determined in 
order to ascertain the presence of absorption and emission features. 
In cases when the dust shell is optically thick at the near-infrared, 
we find that the peak of the pseudocontinuum is at -20 pm, which 
corresponds to a dust temperature of -145 K. The spectral features 
of our sample are compared with pure emission spectra, like those 
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of NGC 6302 and AFGL 4106 [3], in order to determine the shape 
of the absorption features, and the optical depth at different wave- 
lengths. The identification of the several crystalline features is 
done using laboratory spectra. 

Results: We find that the amorphous silicates show extremely 

broad absorption bands; depending on what pseudocontinuum one 
defines, the 18-pm band is well in absorption until approximately 
25-30 jam. Next to amorphous silicates, another material, like 
FeO, may contribute to this absorption feature. On the other hand, 
crystalline silicates occur already in emission at 25 jam and 
longward. Depending on which pseudocontinuum is selected, the 
crystalline silicate features may be in emission, superposed on the 
broad amorphous silicate absorption band. This implies that the 
crystalline and amorphous dust grains may not have the same 
spatial distribution, possibly due to a different radial distribution 
or due to deviations from spherical symmetry. 

References: [1] Waters L. B. F. M. et al. (1996) Astron. 

Astrophys. Lett, 315 , L361. [2] Bedijn P. J. (1987) Astron. Astro - 
phys. y 186, 136. [3] Molster F. J. et al. (1999) Astron. Astrophys ., 
in press. 
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UTILIZING NIGHT SPECTRA OF MARS FOR MINER- 
ALOGY. L. E. Kirkland 1 and K. C. Herr 2 , 1 Lunar and Planetary 
Institute and Rice University, Houston TX, USA (kirk! and @lpi. 
jsc.nasa.gov), 2 The Aerospace Corporation, El Segundo CA, USA. 

Summary: Thermal infrared spectra measured of the night 

side of Mars have advantages over day measurements for exam- 
ining the aerosol dust mineralogy because they contain fewer at- 
mospheric gas and surface features. Despite the advantages night 
spectra offer, they have rarely been utilized because the method- 
ology to convert them to transmission was poorly understood. 
Here we present a method for converting night and day spectra to 
apparent transmission. We also discuss the advantage for miner- 
alogical studies of apparent transmission vs. apparent emissivity. 

Background: In 1969, IRS returned high-quality spectra of 

Mars covering the wavelength region 1.8-14.4 pm. These spectra 
have recently been recovered from the original data tapes and 
calibrated from 5-14.4 pm with heavy involvement from members 
of the original IRS team [1-3]. We use these spectra here because 
IRS has the high signal-to-noise ratio needed for night measure- 
ments. 

Previous Aerosol Studies: Previous studies have used only 

day spectra from IRS or 1971 Mariner IRIS to examine the com- 
position of the aerosol dust. These studies presented a wide vari- 
ety of results, concluding that the the aerosol dust consists mainly 
of an oxidized tholeiitic basalt [4], rock of intermediate [5] or 
felsic composition [6], montmoril Ionite [7], a mixture of basalt and 
clay [8], or palagonite [9]. 

All those studies concluded that the aerosol consists of par- 
ticles <10 pm in size. Widely separated particles this small behave 
in a manner similar to a gas [ 10-12), so their transmission (x) 
and emission (£) spectra have the same shape, although bands that 
appear as transmission troughs appear as emission peaks (x = 
1 - e). 

Night Spectra Advantages: Night measurements contain less 

contribution from atmospheric gases and the surface, which allows 
a more direct examination of the aerosol signature. Night spectra 


have received less attention mainly because it was not clear how 
to convert them from radiance to transmission. 

Since the atmosphere is so cold at night, night spectra contain 
very little contribution from atmospheric hot bands (which are 
from transitions from an excited state rather than the ground state). 
The population of higher energy states decreases exponentially 
with temperature (Boltzmann distribution), so the intensity of a hot 
band decreases rapidly with decreasing temperature. 

Apparent Transmission: Mineralogical studies using thermal 

infrared spectra of Mars typically use spectra converted from ra- 
diance to apparent emissivity (e A ) [10]. However, s A spectra have 
two undesirable artifacts: (1) The spectral contrast depends on the 
background temperature for both day and night spectra, even 
though a true measure of transmission does not have this depen- 
dence; and (2) restrahlen bands in night e A spectra have an emis- 
sivity >1, which is clearly incorrect. 

Apparent emissivity is calculated by dividing the spectrum 
measured of Mars by the Planck radiance curve at the surface 
temperature, taken as the brightness temperature (T B ) at 7.75 pm. 
Figures 1 and 2 show typical day and night IRS spectra, with the 
surface 7.75-jom T B curve shown as a solid line. 

Day measurements are dominated by atmospheric transmission 
and surface emission features, which have restrahlen bands (e.g., 
the 9-pm band) that appear as troughs relative to the surface 
7.75-jom T B Planck curve (Fig. 1). 

However, night spectra have restrahlen bands that appear as 
peaks relative to the surface 7.75-jom T B Planck curve (Fig. 2). 
Thus night e A spectra appear “upside down” relative to day spec- 
tra, and have e A > 1. This is because single-particle reemission 
features dominate night spectra, so restrahlen bands appear as 
peaks rather than troughs. An apparent emissivity calculation fails 
glaringly in this case because it does not account for the reemis- 
sion of the dust. Since the process involved is highly nonlinear, 
the problem cannot be fixed simply by inverting the night e a 
spectra. 

To address this problem, we use a method similar to [13] to 
convert the spectra from radiance to apparent transmission (x A ) 

Mars - aerosol BB 

x A = 

background^ - aerosol BB 

and 

x Aq = 1 - {(1 - x) x log l0 [10 x cos©]} 

where x Aq = x a normalized for atmospheric path length; Mars = 
radiance spectrum measured of Mars; aerosol BB = Planck radiance 
curve at the temperature of the aerosol dust; source BB = Planck 
radiance curve at the temperature of the surface; 0 = emission 
angle. 

The log term in the x Aq equation accounts for absorptions add- 
ing logametrically rather than linearly (Beer’s Law). The surface 
temperature is taken as the brightness temperature at 7.75 jam. The 
aerosol temperature is derived iteratively, and incorporates the 
brightness temperature from spectra measured at the terminator, 
where the flat spectral shape indicates the aerosol and surface 
temperatures are equal. 
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Fig. 1. Conversion of day spectrum to apparent transmission. This shows 
a typical IRS spectrum measured of Mars during the day. To convert to 
apparent transmission, the spectrum is referenced to both the Planck 
blackbody curve of the background (solid line) and the zero reference line 
(dashed line). The background is the surface, which is assumed to behave 
as a perfect blackbody with an emissivity of 1 The zero reference line is 
set at the assumed temperature of the aerosol, and is value that would be 
seen if the transmission of the aerosol is zero. 
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Fig. 2. Conversion of night spectrum to apparent transmission. This 
shows a typical IRS spectrum measured of Mars at night. The spectrum 
is converted to apparent transmission as described above. However, in 
night spectra, the background (surface, solid line) is colder than the zero 
reference line (dashed). The zero reference line gives the value that would 
be seen if the aerosol transmission is zero. 


Discussion: The conversion to x A produces a more accurate 

spectral shape because it accounts for the reemission of the aero- 
sol, thus allowing a more accurate comparison to laboratory spec- 
tra for mineralogical studies of the aerosol. It also provides a 
method to utilize night spectra for mineralogical work. However, 
it assumes the emissivity of the surface is 1, so it does not work 
well where there are strong surface emissivity features. Nonethe- 
less, it works well for Mars, because transmission features domi- 
nate over the very subtle surface features. It is especially suited to 
night measurements, which have even less contribution from the 
surface. 


References: [1] Forney P. B. and Kirkland L. E. (1997) LPS 

XXVIII, 373. [2] Kirkland et al. (1998) LPS XXIX, Abstract #1516. 
[3] Kirkland et al. (1999) LPS XXX, Abstract #1693. [4] Herr 
K. C. et al. (1971) BASS, 3, 466. [5] Hanel R. et al. (1972) Icarus, 
17, 423. [6] Aronsen J. R. and Emslie A. G. (1975) JGR, 80, 
4925. [7] Hunt G. et al. (1973) Icarus, 18, 459. [8] Toon O. B. 
et al. (1977) Icarus, 30, 663. [9] Clancy et al. (1995) JGR, 100, 
5251. [10] Conel J. (1969) JGR, 74, 1614. [11] Hunt G. R. and 
Logan L. M. (1972) Appi Optics, 11, 142. [12] Logan et al. 
(1973) JGR, 78, 4983. [13] Polak M. et al. (1995) Appi Optics, 
34, 5406. 
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MODELING THE INFRARED EMISSION OF DUST IN 
BORN-AGAIN PLANETARY NEBULAE. J Roller 1 and 
S. Kimeswenger 1 , 1 Institute of Astronomy, University of Inns- 
bruck, Technikerstrasse 25/8, A-6020 Innsbruck, Austria (Josef. 
Koller@uibk.ac.at). 

Introduction: We present a model of the IR emission of dust 

in born-again planetary nebulae (PNe) including the emission of 
small dust particles with temperature fluctuations. The IR emission 
in PNe was previously described using a single grain size and a 
temperature equilibrium state (see Pottasch et al. [1] and refer- 
ences therein). This was sufficient to explain the IRAS measure- 
ments. However, this simple equilibrium model cannot explain the 
emission at wavelengths below 10 pm. Bom-again PNe (late He 
flash objects described in Iben et al. [2]) radiate at those wave- 
lengths even stronger than “regular” PNe [3,4). Thus the emission 
has to be described with small dust particles that undergo tempera- 
ture fluctuations. We adopted the technique introduced by 
Guhatha-kurta and Draine [5] with a further detailed description 
by Sieben-morgen et al. [6]. Our results show how important such 
calcula- tions are in the case of a low-density radiation field with 
high-energy photons. Central stars of planetary nebulae (CSPN) 
provide such a dilute field. The required data about the radiation 
field of a late-He flash star, which differs from that of a normal 
CSPN, has been made available by Rauch [7]. 



wavelength [pm] 

Fig. 1. The total emission of carbonaceous dust grains at a certain 
distance from the exciting source. While the two models (solid line = 
transiently heated dust, dashed curve = equilibrium temperature) are very 
similar at longer (IRAS) wavelengths, they differ strongly below 10 pm. 
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Method: The emission of large dust grains (diameter a > 
10 nm = 500,000 C atoms) is calculated with assuming an equi- 
librium between the incident UV radiation and the cooling by 1R 
emission. This can be applied if the cooling timescale x c is sig- 
nificantly larger than the timescale Tj between the incident photons. 
Due to a weak dependency of the absorption cross section, the 
resulting temperature (and the corresponding emission) is only 
weakly dependent on the gram size. A model using a single grain 
size with a single temperature may be applied. 

A small dust grain (a = 0.7 nm = 180 C atoms), which has a 
low heat capacity, is excited by a single high-energy photon to 
rather high temperatures. As the cooling rate is much higher than 
the rate of incident photons, such grains are not in equilibrium. 
The result of a high-energy photon collision is the instant raise of 
temperature, high above the theoretical equilibrium temperature. 
During the cooling process the energy is re-radiated at near infra- 
red wavelengths. The density of the radiation field and its energy 
distribution effects a transition zone from equilibnum temperature 
to transiently heated grains. 

Input Parameters and Results: The spectrum of emitted 

radiation depends on the grain material and the grain size. We use 
carbonaceous grains as in Harrington et al. [3] and the optical 
parameters calculated by Laor and Draine [8]. We show that the 
slope of the grain size distribution (a power law distribution like 
that of Mathis et al. [9]) does not affect the spectrum at wave- 
lengths below 10 |im, while the cutoff radius a min strongly influ- 
ences this emission. Furthermore, we show how the distance from 
the exciting source (thus the density of the radiation field) and the 
spectrum of the CSPN radiation field affects the MIR and FIR 
spectrum. In the end we apply this method to the PN Abell 58 
(= Nova Aql 1919) to fit the unusual IR radiation of this source 
found by ISO [4]. As the different parameters give their “finger- 
prints” to different parts of the emission spectrum, we are able to 
“pin down” the physical properties independently. 
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ISOPHOT POLARIZATION MEASUREMENTS OF THE 
ASTEROIDS 6 HEBE AND 9 METIS AT 25 MICRO- 
METERS. J. S. V. Lagerros 1 , T. G. Muller 2 , U. Klaas 3 , and A. 
Erikson 4 , ‘Five College Radio Astronomy Observatory, LGRC- 
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Fig, 1. Model thermal lightcurves for 6 Hebe at 25 pm, assuming a 
refractive index of n = 2, and r.m.s. slopes p = 0, 0.1, . . . , 1.0. 
(a) Observed and predicted flux / v . The model fluxes increase with the 
surface roughness p. To the right is the typical error bar for the model 
fluxes, due to the uncertainty in the occultation diameter only, (b) The 
time span and detection degree of linear polarization 7t L . The predicted 
degree of linear polarization increases with lower surface roughness. The 
upper limit by the ISOPHOT observation is shown. 


The subsurface thermal IR emission of asteroids becomes lin- 
early polarized when scattered across the surface. The disk inte- 
grated mid-IR Stokes parameters are computed by extending the 
new thermophysical model by Lagerros [1,2]. The small hemi- 
spherical segment craters covering a smooth surface. Analytical 
solutions are used for the multiple scattered solar and thermally 
emitted radiation inside the craters. The surface roughness en- 
hances the emission toward the solar direction, but randomizes the 
scattering planes which lowers the polarization. In general the 
predicted degree of linear polarization increases with higher refrac- 
tive index of the surface material, and a more elongated shape of 
the asteroid. 

The asteroids 6 Hebe and 9 Metis were observed at X = 25 pm 
with ISOPHOT onboard the Infrared Space Observatory (ISO). 
The model absolute fluxes are in good agreement with the photo- 
metric results. Although no linear polarization was detected, the 
upper limits together with the extended model enable us to put 
useful constraints on the regolith properties of the target asteroids. 
The derived detection limits are compared to model polarization, 
by spanning a range in surface roughness, refractive index, and 
thermal inertia parameter space. The Metis observations favor a 
low refractive index and high surface roughness, but the Hebe 
observations are inconclusive since they coincided with a mini- 
mum in the polarization curve. 

References: [1] Lagerros J. S. V. (1996) Astron. Astrophys ., 
315, 625-632. [2] Lagerros J. S. V. (1998) Astron. Astrophys., 
332, 1123-1132. 



THE DUST EMISSION FEATURES IN THE INFRARED 
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Laboratory, Chilton, Didcot, Oxon 0X11 OQX, UK (T.Lim@ 
rl.ac.uk), 2 Astronomical Institute Anton Pannekoek, University of 
Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 
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Introduction: Spectra in the 2.5-97-pm region obtained by 

the Infrared Space Observatory (ISO) SWS and LWS spectrom- 
eters have provided evidence for the presence of crystalline sili- 
cates in a number of different objects. The first evidence was 
provided by the ISO SWS instrument and was reported soon af- 
ter the ISO launch by Waters et al. [1]. The reason that laboratory 
studies had concentrated on glassy/amorphous silicates prior to 
ISO was because the broad 10-pm silicate features observed from 
the ground and with the IRAS LRS were broad and smooth and 
could not be fitted with crystalline silicates. Since ISO’s discov- 
ery of large concentrations of astronomical crystalline silicates, 
efforts have been made to characterize their optical properties at 
infrared wavelengths, most recently by the group at the University 
Observatory in Jena [2], where new optical constants for the as- 
tronomically important minerals, e.g., enstatite, have recently be- 
come available. 

We report on the study of a range of broad features exhibited 
in the spectrum of the planetary nebula NGC 6302. Full-range 
grating spectra were taken by both ISO spectrometers (SWS and 
LWS) with a combined wavelength coverage 2.5-197 pm. The 
spectrum, one of the richest observed by ISO, exhibits prominent 
broad features between 3.3 and 1 17 pm that are attributed to vari- 
ous silicate minerals, water ice, and PAH modes. A spline fit was 
made to a set of points assumed to be indicative of the continuum 
level and the spectrum was divided by the fit. The species present 
were then identified from the continuum divided spectrum and 
these were characterized by multiple gaussian fitting. An attempt 
has also been made to fit the whole spectrum with a combina- 
tion of silicate species and water ice. The overall fit has two 
problems: The peak wavelengths of the emission bands of some 
species given by the model, which are based on laboratory spec- 
tra of minerals, do not match those observed with ISO, also not 
all the features are fitted. The mismatch of wavelengths can be at- 
tributed to probable differences between the composition and/or 
particle size of the samples used to obtain the optical constants and 
those present in the environment of NGC 6302. The origin of the 
observed features that lack laboratory matches is unknown and 
until optical constants from pure minerals are available in the ISO 
wavelength region, it is not possible to determine whether these 
unidentified features are due to impure samples of the minerals 
used for our fits, or to minerals whose optical constants are not yet 
available. 

References: [1] Waters L. B. F. M. et al. (1996) Astron. 
Astrophys ., 315, L361. [2] Jager C. et al. (1998) Astron. Astro- 
phys., 339 , 904. 
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ADVENTURES IN MODELING THERMAL EMISSION 
SPECTRA OF COMETS. D. K. Lynch and S. Mazuk, The 
Aerospace Corporation, M2/266, P.O. Box 92957, Los Angeles 
CA 90009, USA. 

We review some of the current problems in the interpretation 
of 8-13-jim thermal emission spectra and discuss a number of 


aspects of small particles that are seldom included in modem 
work. Topics include relating observed color temperature to kinetic 
temperature, the effects of size and optical constants on tempera- 
ture, the nearly ubiquitous existence of molecular disorder, estimat- 
ing the optical constants of disordered material, the effects of 
nonspherical particles, and the wavelength shift of lattice reso- 
nances with particle size. We develop the analyses around the 
silicate emission feature observed in many comets and discuss why 
comets are ideal subjects for advancing our knowledge of thermal 
emission from small particles. The results are broadly applicable 
to all cosmic dust. 



INFRARED SPACE OBSERVATORY SPECTROSCOPY 
OF CIRCUMSTELLAR DUST DISKS SURROUNDING 
“ISO”LATED HERBIG Ae/Be STARS, THEIR EVOLU- 
TION, AND THE RELATION WITH OUR OWN SOLAR 
SYSTEM. K. Malfait 1 and C. Waelkens 1 , Hnstituut voor 
Sterrenkunde, Celestijnenlaan 200B, B-3001 Heverlee, 
Belgium (Koen.Malfait@ster.kuleuven.ac.be; Christoffel.Waelkens@ 
ster.kuleuven.ac.be). 

Introduction: During the last decade, there has been a grow- 

ing interest in circumstellar dust disks surrounding young stellar 
objects. Herbig Ae/Be stars are assumed to be the evolutionary 
precursors of Beta Pictoris like stars. Today, the most studied 
Herbig Ae/Be stars are still embedded in the dust envelope from 
which they are formed. The isolated Herbig Ae/Be stars we study 
are not obscured by filaments of a star formation region; all the 
infrared excess they possess comes from their circumstellar dust 
disk. Therefore, these stars are very good targets for infrared spec- 
troscopy. 

We will present the ISO spectra of several isolated Herbig Ad 
Be stars and we will demonstrate that the spectra evolve during the 
pre-main-sequence phase. Remarkable similarities with solar sys- 
tem dust occurs. The ISO spectrum of Beta Pictoris will be men- 
tioned as well. 

Observations and Reduction: In the guaranteed time pro- 

posal of Waelkens, a sample of more than 20 isolated Herbig Ad 
Be stars with reasonable infrared fluxes has been studied using the 
Short Wavelength Spectrometer onboard the Infrared Space Obser- 
vatory (ISO). In five cases, this 2.3-45-pm spectrum is combined 
with a spectrum obtained with the Long Wavelength Spectrometer. 
The resulting resolution of the spectra varies from 250 to 1500, 
depending on the target dedicated time. 

Results: While the youngest, most embedded Herbig Ae/Be 

stars [1] exhibit ice absorption bands and amorphous silicates, the 
spectra of the isolated Herbig Ae/Be stars seem to be dominated 
by optically thin solid state emission features. A wide variety of 
solid state features is visible, such as the PAH bands, silicates 
(both amorphous and crystalline), oxides, and cold water ice. 

The spectra of some isolated, and therefore assumed to be the 
most evolved, Herbig Ae/Be stars present astonishing similarities 
with solar system objects. The ISO spectrum of HD 100546 re- 
veals a forest of solid-state features attributed to crystalline sili- 
cates, mainly forsterite, the Mg-rich end member of the olivine 
family. Next to this crystalline dust component, also amorphous 
silicates, FeO, PAHs, and crystalline water ice are present [2]. The 
similarity between the spectrum of HD 100546 and that of Comet 
Hale-Bopp [2-4] is astonishing The crystalline material in the 
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environment of HD 100546 and in the comet has a temperature of 
less than 250 K. However, laboratory studies suggest that silicates 
can only crystallize at temperatures of more than 1000 K [e.g., 5]. 

The hot passages close to the Sun are too short to be able to 
crystallize the Hale-Bopp material. Therefore, we argue that the 
dust must have crystallized closer to the hot stellar photosphere. 
This implies that the radiation we see in the environment of HD 
100546 emerges from a forming Kuiper Belt or an Oort Cloud. 
Indirectly, the presence of a giant planet in the environment of HD 
100546 is likely, since the Oort Cloud and a Kuiper Belt in our 
solar system originate from ejection of inner solar system material 
by Neptune and Uranus [6]. 

HD 142527, a star much younger than HD 100546, does not 
exhibit the cold crystalline silicates, although some evidence for 
hot crystalline silicates is present. The cold dust radiation is com- 
posed of three components: Next to the omnipresent amorphous 
silicates and very cold crystalline water ice, hydrous silicates are 
a major constituent of the circumstellar dust [7]. A large fraction 
of the interplanetary dust particles that are gathered during atmo- 
spheric flights in our own atmosphere also have this lattice layer 
structure [8]. The temperature of the crystalline ice and the hy- 
drous silicates again does not correspond to what can be expected 
from laboratory experiments. The location of the hydrous silicates 
further outside the stellar disk than the crystalline silicates is con- 
sistent with their temperature. 

We will also present the spectra of HD 45677 and HD 179218, 
which are remarkably alike, although the nature and evolutionary 
status of HD 45677 is not clear. The age of HD 179218 and HD 
142527 is comparable, though their ISO spectra differ strongly [9]. 

Conclusions: The similar nature of the solar system dust and 

the circumstellar dust surrounding young stars confirms that solar 
system formation occurs at this stage of stellar evolution. Although 
the resemblance between some spectra is very strong, a wide 
variety of spectra occurs. The estimated ages for HD 142527 and 
HD 179218 are of the same order, but their spectra do not re- 
semble each other. 

Although there is no evident correlation between stellar param- 
eters, stellar age and the observed spectra, a general trend can be 
distinguished, going from optically thick radiation, characterized 
by ice absorption features and amorphous silicates, toward opti- 
cally thin radiation that partially originates from crystalline dust. 
The 10-pm silicate band of Beta Pictoris is dominated by crystal- 
line silicates as well [10,11]. This confirms that the isolated Herbig 
Ae/Be stars are in an evolutionary passage from their embedded 
analogs toward Beta Pictoris objects, as can be seen from the 
evolution of the spectral energy distributions as well [12]. 

References: [1] Waelkens C. et al. (1998) to be published in 

the proceedings of the “First International Conference on Comet 
Hale-Bopp,” held in Tenerife, March 1998. [2] Malfait K. et al. 
(1998) Astron. Astrophys., 331, L25— L28. [3] Crovisier J. et al. 
(1997) Science, 1997 , 1904-1907. [4] Lellouch E. et al. (1998) 
Astron. Astrophys ., 339 , L9-L12. [5] Hallenbeck S. et al. (1998) 
Icarus, 131, 198-209. [6] Fernandez J. A. (1978) Icarus, 34, 173- 
181. [7] Malfait K. et al. (1999) Astron. Astrophys., submitted. 
[8] Sandford S. A. and Walker R. M. (1985) Astrophys. J., 291, 
838-851. [9] Waelkens C. et al. (1998) to be published in the 
proceedings of “The Universe as seen by ISO,” a conference held 
in Paris, France, October 1998. [10] Knacke R. F. et al. (1993) 
Astrophys. J., 418, 440—450. [11] Pantin C. K. (1998) to be pub- 
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A COMPARATIVE STUDY OF INFRARED SPECTRAL 
FEATURES IN HERBIG Ae/Be STARS. G Meeus 1 and C 
Waelkens 1 , Astronomical Institute, Celestijnenlaan 200B, B-3001 
Heverlee (gwendolyn@ster.kuleuven.ac.be). 

Introduction: The observations of the Infrared Space Obser- 

vatory (ISO) satellite revealed an amazing variety of spectral fea- 
tures in young stars. Both O- and C-rich dust was found, and 
silicates were seen both in amorphous and crystalline form. An 
important conclusion from ISO observations is that these crystaline 
silicates are, apart from seen in young stars [1], also found in 
meteors, comets [2], evolved stars, and interplanetary dust par- 
ticles. We studied ISO spectra for a sample of 30 young stars, of 
which most are Herbig Ae/Be stars. Correlation between dust 
features on the one hand and extinction, age, spectral type, and 
amount of IR excess on the other are looked for. We notice an 
evolution in composition and amount of the dust when the star 
gets closer to the Main Sequence. 

Observations: We dispose over ISO Short Wavelength Spec- 

trometer (SWS) observations of a large sample of isolated Herbig 
Ae/Be stars, covering a wide range in spectral type. These spec- 
tra are part of the guaranteed time proposals of Waelkens and 
Wesselius. Observations are done with ISO SWS in mode AOT1; 
the spectra cover an interval from 2.3 to 45 pm. For five of the 
stars we also have LWS spectra. The resolution varies between 
250 and 1500, depending on the integration time. 

Results: Classification scheme. The most important star-to- 

star differences are the absence/presence of silicate features, cool 
dust and PAHs, and a rising/flat continuum. From these differ- 
ences seen in the spectra, we can distinguish five main groups: (1) 
rising, featureless spectrum (embedded sources); (2) rising, rela- 
tively featureless, silicate in absorption or absent; (3) rising with 
a prominent cool dust component, silicate emission; (4) silicate 
emission without a prominent cool dust component; and (5) fea- 
tureless without cool dust. Further subdivision can be made, de- 
pending on the presence of PAHs and on the structure of the 
silicates. 

Comparison with optical observations. We have already stud- 
ied most of our sources through optical spectra and photometry. 
This studies reveal some of the sample stars to be variable, have 
outflows, or show signs of metal underabundance. 

These findings were confronted with properties of the CS dust. 
Correlation between strength, structure, and presence of dust fea- 
tures on the one hand and spectral type, E[B-V], presence of ac- 
cretion/outflow, chemical composition and evolutionary stadium on 
the other were looked for. 

The following results were found: (1) There is no relation be- 
tween the strength of the silicate features and the extinction (E[B- 
V]). (2) There is no relation between the spectral type of the star 
and presence/absence, structure, or strength of spectral features. 
(3) All stars showing crystalline silicates have a prominent cool 
dust continuum, and also show PAHs. PAHs are only present in 
stars from group 1 to 3. (4) Disk inclination has no influence on 
the IR spectral features, suggesting the material to be optically 
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thin. (5) No activity (infall or outflow) is seen in group 5. (6) Stars 
that are underabundant in metal all belong to group 5, and there- 
fore are featureless. The stars with a solar composition show sili- 
cates in emission. 

Conclusions: Group 1 stars have featureless IR spectra and 

are embedded, thus they are the youngest of the whole sample. 
Group 5 stars are inactive and not embedded, thus they probably 
are the most evolved ones. 

These results suggest that the order of the groups indicates an 
evolutionary sequence and that PAHs and silicates seem to disap- 
pear with age. 

The absence of PAHs and silicates can be explained by a lower 
UV flux (in the case of the PAHs), and by the formation of larger 
bodies out of small particles (in the case of the silicates). 

As a result, the oldest stars from our sample are suggested to 
be the closest progenitors of B Pic and Vega-type stars. 

References: [1] Malfait K. et al. (1998) Astron. Astrophys 
331 , L25-L28. [2] Crovisier J. et al. (1998) Science , 1997, 1904- 
1907. 



CRYSTALLINE SILICATES: FROM NOTHING TO MAIN 
THING. F. J. Molster 1 and L. B. F. M. Waters 1 - 2 , Astronomical 
Institute Anton Pannekoek, University of Amsterdam, Kruislaan 
403, 1098 SJ Amsterdam, The Netherlands, 2 SRON Space Re- 
search Laboratory, P.O. Box 800, 9700 AV Groningen, The 
Netherlands (frankm@astro.uva.nl). 

Introduction: One of the biggest surprises of the Infrared 

Space Observatory (ISO) was the discovery of crystalline silicates 
in the outflows of O-rich evolved stars. Before the launch of ISO 
it was commonly accepted that the silicates, apart from a few 
exceptions in the remnants of some accretion disks of young stars 
and of our own solar system accretion disk in the form of com- 
ets and interplanetary dust particles (IDPs), were of amorphous 
structure. The spectra taken with the Short and Long Wavelength 
Spectrometer (SWS/LWS) onboard ISO showed clear evidence of 
the presence of a significant fraction (5-15%) of crystalline sili- 
cates in the circumstellar environment (CSE) of evolved stars [1]. 
The most common found crystalline silicates are the (almost) pure 
Mg end members of the olivine (Mg 2x Fe (2 _ 2x) Si0 4 ) and pyroxene 
(Mg x Fe {1 _ x) Si0 3 ) solid solution senes, i.e., x > 0.99, and are called 
forsterite (Mg 2x Si0 4 ) and enstatite (Mg x Si0 3 ) [1,2]. This is in 
contrast with the relatively Fe-rich amorphous silicates. Since the 
Mg-silicates condense at higher temperatures than the Fe-rich sili- 
cates (see, e.g., Gail and Sedlmayr, in press) it is expected that the 
crystalline silicates are formed in the very early stages of the dust 
formation. Differences in the appearance of the crystalline silicates 
are therefore probably related to the different conditions in the dust 
forming areas. 

Discussion of the Observations: In order to understand more 

of the (crystalline) dust formation we have observed with ISO 
SWS and LWS several evolved stars with crystalline silicate fea- 
tures. Where most stars have modest crystalline silicate abundances 
(<15%), some stars show a significantly higher abundances, up to 
75%, of crystalline silicates. It appeared that all stars with a high 
abundances of crystalline silicates have a disk. Apart from the fact 
that the abundance of crystalline silicates might be a new indica- 
tor for the presence of a disk, it is also very intriguing how these 


disks become so highly crystalline. The relatively low abundance 
found in “normal” outflow sources suggests that extra processing 
must have taken place in the disk. This is not an obvious process, 
since not only should most of the dust particles become crystal- 
line, the chemical structure of the particles also has to change from 
Fe-rich to Fe-free. We shall discuss different scenarios to produce 
such highly crystalline dust disks in the context of these and other 
ISO observations. 

References: [1] Waters L. B. F. M. et al. (1996) Astron. 

Astrophys., 315, L361. [2] Jager C. et al. (1998) Astron. Astro- 
phys., 339, 904. 
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MARTIAN MINERALOGY WITH THE INFRARED SPACE 
OBSERVATORY SHORT WAVELENGTH SPECTRO- 
METER. P. W. Morris 1 ’ 710 , T. de Graauw 2 , E. Lellouch 3 , B. 
Henderson 4 , T. Encrenaz 3 , S. Erard 5 , H. Feuchtgruber 6 , M. 
Burgdorf 7 , G. Davis 8 , and M. Griffin 9 , l SRON Laboratory for 
Space Research, Sorbonnelaan 2, NL-3584 CA Utrecht, The 
Netherlands, 2 SRON Laboratory for Space Research. Landlevan 
12, NL-9747 AD Groningen, The Netherlands, 3 DESPA, 
Observatoire de Paris, Meudon, France, 4 Earth Watch Inc., 1900 
Pike Road, Longmont CO 80501, USA. 5 CNRS/Umversite Paris- 
Sud bat. 121, F-91405 Orsay, France, 6 MPE Garchmg. D-85748 
Garching, Germany, 7 ISO Science Operations Center, Vilspa, E- 
28080 Madrid, Spain, 8 University of Saskatchewan, Saskatoon, 
Saskatchewan, Canada, 9 Queen Mary and Westfield College, 
London El 4NS, UK, 10 Astronomical Institute, University of 
Amsterdam, NL-1098 SJ Amsterdam, The Netherlands. 


Thermal infrared spectroscopic observations of Mars covering 
2.4-45 pm were obtained in mid- 1997 with the Short Wavelength 
Spectrometer (SWS) onboard ESA’s Infrared Space Observatory. 
At very high signal-to-noise ratios and spectral resolutions of up 
to 3000 using the grating sections of SWS, these data provide 
composite views of the gas and solid-state chemistry in three dif- 
ferent pointings with unprecedented detail. Three observations with 
complete wavelength coverage were obtained in July— August 
1997. with subsolar points located close to Elysium-Isidis Plantia, 
Tharsis-Lunae Planum, and Syrtis Major. We present our progress 
in analyzing these observations, with special emphasis on the 
features attributed to minerals of the dust aerosol and surface. A 
synthetic spectrum consisting of the surface thermal emission and 
gas atmosphere has been synthesized for each observation, using 
climate information from the ESA Martian Climate Database and 
the observed C0 2 and H 2 0 vapor lines. Division of the observed 
spectra by the synthetic atmospheres reveals an array of features 
that have allowed us to make certain mineral identifications by 
comparison to available libraries of laboratory mineral reflectance 
measurements. In this study we address the particulate effects of 
near-surface thermal gradients and resulting influence on band 
positions with respect to isothermal conditions and bulk samples. 
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ASTEROID 4 VESTA AS SEEN WITH THE INFRARED 
SPACE OBSERVATORY SHORT WAVELENGTH SPEC- 
TROMETER. P. W. Morris 1 - 2 , A. M. Heras 3 , B. Vanden- 
bussche 4 , and R. Dijkstra 1 , Astronomical Institute, University of 
Amsterdam, The Netherlands (pmorris@astro.uva.nl), 2 SRON 
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Laboratory for Space Research, Utrecht, The Netherlands, Euro- 
pean Space Agency, ESTEC/SCI-SAF, Noordwijk, The Nether- 
lands, institute of Astronomy, Catholic University of Leuven, 
Belgium. 

The geological uniqueness and diversity of asteroid 4 Vesta is 
known from rotationally resolved spectroscopic measurements in 
the visible and near- infrared, and from polarimetry measurements. 
These measurements have demonstrated longitudinal variations in 
surface composition and albedo that can account for its lightcurve, 
compared to normal asteroid hghtcurves induced by shadowing. 
Olivine in at least one small region has been previously detected, 
and is thought to have been excavated from the ferromagnesian 
mantle by a major impact event. The impact region is obvious in 
1994 HST images. On the basis of the hemispheric dichotomy of 
albedo apparent in the HST WFPC2 images, the impact would 
seem not to have resurfaced the entire asteroid. 

We have obtained a series of high spectral resolution (1000- 
3000) scans of 4 Vesta over the 4-45-pm range with the ISO 
Short Wavelength Spectrometer. The scans are resolved in steps 
of one quarter of Vesta’s rotation period, of -24-min duration 
(0.075 P). In this presentation, we report the detection and vari- 
ability of strong Si-O-Si asymmetric stretch fundamental of oliv- 
ine near 11 pm (emission), and tentative pyroxene near 10 pm 
(absorption). 

EFFECTS OF GRAIN MORPHOLOGY AND IMPURITIES 
ON THE INFRARED SPECTRA OF SILICON CARBIDE 
PARTICLES. H. Mutschke, D. Clement, and Th. Henning, 
Astrophysical Institute and University Observatory, Schiller- 
gasschen 2-3, D-07745 Jena, Germany (mutschke@astro.uni- 
jena.de). 

Introduction: Thermal emission of solid silicon carbide (SiC) 

particles has been observed in carbon star outflows [1-3]. The 
identification is based on a broad infrared feature in the 10-1 3-pm 
wavelength range, which is ascribed to the fundamental phonon 
transition of SiC. Attempts have been made to derive the crystal 
type of the circumstellar SiC grains from a comparison of the 
observed band profiles to those measured on laboratory analogs 
[3]. However, due to the strength of the phonon transition the 
absorption by small SiC grains is determined by surface modes, 
which make the band profile extremely sensitive to the grain 
morphology (size, shape, and grain aggregation [4]). Unfortu- 
nately, at the moment light scattering theory is not in the position 
to calculate band profiles for arbitrary grain morphologies. There- 
fore, intensive laboratory studies must help to clarify the influence 
of morphology on the appearance of the SiC phonon band. 

Experimental: We present experimental infrared absorption 

spectra, measured on SiC particulate samples of different grain 
sizes and shapes as well as in different agglomeration states. For 
this purpose, a number of preparational methods have been ap- 
plied, including isolation in noble gas matrices, the KBr pellet 
technique and IR microscopy of grains supported by Si substrates. 
For the matrix-isolation spectroscopy, SiC particles have been 
produced in a gas-pyrolysis source. 

Results and Discussion: With SiC particulate samples of 

different (partly commercial) ongin we observe a wide variety of 



Fig. 1 . Phonon band profiles of some SiC laboratory samples. The 
spectra have been shifted vertically for clarity, a and p indicate the 
hexagonal/rhombohedral and the cubic crystal types respectively. 


SiC phonon band profiles, which is not primarily related to the 
crystal type (see Fig. 1). We show that both peak wavelength(s) 
and band shape depend on grain size, shape, and agglomeration 
state, as well as on the medium surrounding the grains. For a more 
detailed analysis of these effects it is necessary to separate the 
agglomeration effect from the one produced by the primary grain 
morphology. Matrix isolation can provide a means to solve this 
problem in the future. 

Additionally, we find that surface plasmon absorption caused 
by impurities influences the phonon band via plasmon-phonon 
coupling. Moreover, it produces a very strong quasicontinuous 
absorption in the near- and mid-infrared, which could be impor- 
tant to take into account in models of the thermal dust emission 
of carbon star envelopes. 

References: [1] Little-Marenin 1. R. (1986) Astrophys. J., 

307 , L15. [2] Goebel J. H. et al. (1995) Astrophys. J 449, 246. 
[3] Speck A. K. et al. (1997) MNRAS, 288, 431. [4] Papoular R. 
(1998) Astron. Astrophys., 329, 1035. 
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DUSTY: A PUBLICLY AVAILABLE CODE FOR CON- 
TINUUM RADIATIVE TRANSFER IN ASTROPHYSICAL 
ENVIRONMENTS. M. Nenkova 1 , Z. Ivezic 2 , and M. Elitzur 1 , 
apartment of Physics and Astronomy, University of Kentucky, 
Lexington KY 40506-0055, USA (maia@pa.uky.edu; moshe@ 
pa.uky.edu), 2 Astrophysical Sciences Department, Payton Hall, 
Princeton University, Princeton, NJ 08544-1001, USA (ivezic @ 
astro.princeton.edu). 
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External radiation . Thanks to scale invariance [2], the radia- 
tion illuminating the dusty shell needs to be specified only by its 
spectral shape. DUSTY provides options for entry in analytical or 
numerical form, as user-supplied files. 

Dust optical properties. These are described by the absorp- 
tion and scattering cross sections, which depend on the grain size 
and material. DUSTY contains data for the optical properties of 
six common grain types: Ossenkopf ‘‘warm” and “cold ’ silicates 
[4], silicates and graphite grains from Draine and Lee [5], amor- 
phous carbon from Hanner [6], and SiC from Pegourie [7]. Ad- 
ditional optical properties for other grains can be supplied if 
needed. In this case, the user can either specify directly the absorp- 
tion and scattering coefficients or have DUSTY calculate them 
from provided optical constants by using the Mie theory. Two 
types of size distribution functions are implemented in the 
code: MRN [8] and its modification KMH [9]. The user is free 
to adjust the parameters of size distribution as desired. 

Density distribution . DUSTY accepts three methods for en- 
tering the density distribution: prescribed analytical forms, hydro- 
dynamic calculation of winds driven by radiation pressure on dust 
particles, and numerical tabulation in a file. 

Optical depth . For a given set of the parameters specified 
above, DUSTY can generate up to 99 models covering a range of 
optical depths x 0 , specified at some fiducial wavelength Kq. 

Output. A typical DUSTY run generates an enormous amount 
of information, and the user is in control of the volume of output, 
which may include spectra, imaging profiles, and radial distribu- 
tions for each of the optical depths included in the run. 

Details and instructions about the many options can be found 
in DUSTY’s manual [1]. The code is publicly available for vari- 
ous modeling purposes and the authors would appreciate any 
comments and suggestions leading toward improvement of the 
code and its capabilities. 

References: [1] Ivezic Z. et al. (1999) User Manual for 

DUSTY, Internal Report, Univ. of Kentucky (accessible at http:// 
www.pa.uky.edu/~moshe/dusty). [2] Ivezic Z. and Elitzur M. 
(1997) MNRAS, 287 , 799. [3] Schmidt-Burgk J. (1975) Astron. 
Astrophys., 40, 249. [4] Ossenkopf V.et al. (1992) Astron. 
Astrophys ., 261, 567. [5] Draine B. and Lee H. (1984) Astrophys. 
J., 285, 89. [6] Hanner M. (1988) NASA Conf. Pub. 3004, 22. 
[7] Pegourie B. (1988) Astron. Astrophys 194, 335. [8] Mathis 
J. et al. (1977) Astrophys. J., 217, 425. [9] Kim S. et al. (1994) 
Astrophys. J., 422, 164. 



NUCLEATION, GROWTH, ANNEALING, AND COAGU- 
LATION OF REFRACTORY OXIDES AND METALS: 
RECENT EXPERIMENTAL PROGRESS AND APPLICA- 
TIONS TO ASTROPHYSICAL SYSTEMS. J A Nuth 1 , 
F. J. M. Rietmeijer 2 , S. L. Hallenbeck 3 , and P. A. Withey 4 , 
l Astrochemistry Branch, Mail Code 691, NASA Goddard Space 
Flight Center, Greenbelt MD 20771, USA, 2 Geology Department, 
University of New Mexico, Albuquerque NM 85131, USA, 3 Mail 
Code SR, NASA Headquarters, Washington DC 20546, USA, 
department of Physics and Engineering, West Virginia Wesleyan 
College, Buckhannon WV 26201, USA. 

Introduction: Starting with cooling, refractory vapors diluted 

in significant quantities of H and He there are four processes that 


most natural systems will undergo: nucleation, growth, annealing, 
and coagulation. Nucleation is the processes by which the first 
stable refractory nuclei form in the vapor. These are the seeds onto 
which the remaining vapors will condense during the growth stage. 
Solids of any composition will try to arrange themselves into the 
least energetic configuration, provided that there is sufficient en- 
ergy available to support such processes as diffusion and the 
breaking of chemical bonds. There is a significant activation en- 
ergy associated with the annealing process in refractory solids due 
to the relatively high energy of the chemical bonds in solids. The 
grains formed in most cosmochemical systems are extremely small 
and often tightly coupled to the gas. Because of their small physi- 
cal cross sections coagulation may be a very slow process unless 
there is another driving force involved in addition to normal 
Brownian motion. In what follows we will briefly cover each of 
these four stages for refractory oxide and metal grams, although 
in inverse order. 

Coagulation: Appropriate studies of the coagulation of re- 

fractory grains are just beginning to be carried out. In general, 
such studies are inherently difficult since, just as the grain size 
begins to become interesting, the aggregate settles out of the gas. 
To overcome this problem, experimenters are beginning to take ad- 
vantage of the short-lived periods of microgravity available in drop 
towers and aircraft in preparation for work to be carried out on the 
International Space Station. We have recently demonstrated that 
the nucleation of very small Fe grains can lead to the formation 
of single-magnetic domain particles, e.g., fully magnetized iron 
dipoles [1] that might greatly enhance the coagulation efficiency 
of small particles [2]. This could lead to the formation of very low 
fractal dimension particles that might greatly increase the emission 
efficiency of these aggregates at submillimeter wavelengths [3). 

Annealing: Over the last several years we have made signifi- 

cant progress in understanding spectra of circumstellar and inter- 
stellar oxide grains. In particular, based on experimental studies in 
our lab [4] we can now predict the mid-IR spectral evolution of 
silicate grains as they anneal in stellar outflows and in forming 
protostellar nebulae. This not only allows us to predict the emer- 
gent spectrum one might expect from obscured sources, but also 
allows us to set limits on the time-temperature history of dust 
grains based on their spectral properties. As an example, we have 
been able to show that the spectral properties of grains in Hale- 
Bopp, Halley, and other comets are the result of moderate tem- 
perature (1000 K) thermal processes in the solar nebula that 
require recirculation from the inner nebula to well beyond Jupiter’s 
orbit [5]. These large-scale circulation patterns are not currently 
part of any model of protostellar evolution, but, if they exist for 
any significant timescale, would have a large effect on the nebu- 
lar production efficiency of prebiotic organic compounds and on 
the composition of volatile compounds sequestered in comets. 

Growth: Another area where we have made significant pro- 

gress in the last two years is in understanding the chemical spe- 
ciation of grains as they grow in stellar outflows. Data acquired 
by careful analysis of the chemical compositions of individual 
grains condensed from mixed refractory vapors in our lab has 
demonstrated that these compositions are predictable based on the 
compositions of metastable eutectics in the appropriate phase dia- 
grams [6—8] such as that shown in Fig. 1. Based on these experi- 
ments, we now hypothesize that the chemical compositions of the 
individual grains formed in such outflows are controlled by ther- 
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Introduction: Many astronomical objects are surrounded by 

dust and often can be entirely obscured. The dust scatters, absorbs, 
and reemits the radiation onginating from the source, which results 
in an overall shift of the spectral energy distribution toward longer 
IR wavelengths. In many cases the emerging processed spectrum 
provides the only available information about the dust-embedded 
object. Thanks to the recent developments in IR techniques and 
facilities, the IR signature of many objects is now available, and 
the amount and quality of data is steadily increasing. Interpreta- 
tion of these data can offer insight to the nature of optically ob- 
scured objects, making their IR signature a powerful tool of 
analysis. However, because of the complexity of the radiative 
transfer problem, proper computational modeling tools must aid 
such analysis. We have developed a publicly available code 
DUSTY, which solves the radiative transfer problem for a source 
embedded in a spherically symmetric dusty envelope. The code 
offers also an optional calculation for a dusty plane-parallel slab, 
illuminated from one or both sides at an arbitrary angle. The user 
specifies the properties of the radiation source and dusty medium, 
and the code calculates the dust temperature distribution, radiation 
field, and other quantities of interest [1], 

Solution Method: The code fully utilizes the scaling proper- 

ties of the radiative transfer problem, which are thoroughly dis- 
cussed in [2]. DUSTY solves the integral equation for the energy 
density in a dusty shell obtained by analytically integrating the 
radiative transfer equation. This equation, coupled with the con- 
dition of radiative balance at any point of the envelope, is solved 
by a direct method, i.e., the numerical integration over angles and 
optical depth is transformed into multiplication with a matrix of 
weight factors [3], The matrix is determined purely by the geom- 


etry and is calculated only once for given space grids. The energy 
density and the next radiative moment — the flux at every point 
are then determined by direct matrix inversion, which solves the 
scattering problem and eliminates the need of iterations over the 
energy density itself. That gives our method an edge in speed over 
other schemes (e.g., iterations over variable Eddington factor). The 
employed direct method of solution provides fast and stable con- 
vergence even for optical depths as large as t v -1000. The radia- 
tive transfer problem in spherical symmetry is solved on grids for 
dimensionless radius and impact parameter, whereas for the plane- 
parallel case the calculation is performed in optical depth space. 
The wavelength grid is predetermined in the code and contains 98 
points from 0.01 jam to 3.6 mm. The spatial grids are automati- 
cally generated and adjusted to ensure flux conservation. 

Code Capabilities: One of the goals while developing 

DUSTY was to design a flexible, user-friendly input/output inter- 
face. Another advantage of the code is that the number of inde- 
pendent input model parameters is minimal since it fully utilizes 
the scaling properties of the radiative transfer. The input contains 
three types of data: physical parameters, numerical accuracy pa- 
rameters, and flags for optional output files. The physical param- 
eters include characteristics of the external radiation, properties of 
the dust grains, and the envelope density distribution. An example 
of a code application is shown in Fig. 1., a set of models for an 
envelope with silicate dust and 1/r 2 density distribution around a 
2500 K central star. 

Plane-parallel slab. This option offers calculation of the ra- 
diative transfer through a dusty homogeneous slab of material il- 
luminated from one or both sides at user-specified angles of 
incidence (Fig. 2). 


Model Spectra for Density Distribution 
100% Silicates, MRN Size Distribution 
Total Optical Depth Ty= 0.1-1000 




Model Results for a Slab of ISM Dust with Total Optical Depth Ty 
Normal Illumination by a 10,000K Black Body on the Left 
and a 3.000K Black Body on the Right 




Fig. 1. 


Fig. 2. 
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Fig. 1 . Phase diagram for the MgO-Si0 2 system showing the com- 
positions of the observed metastable eutectics that also correspond to 
observed abundance peaks in the composition of the grains condensed 
from mixed Mg-SiO-H vapors. 


modynamic factors and can be predicted by appropriate chemical 
phase diagrams. Although the relative proportions of the individual 
grain compositions will vary from star to star based on the over- 
all chemical composition of the outflow; the compositions of the 
individual grains themselves will remain unchanged. This obser- 
vation explains the remarkable similarity of the spectra of grains 
produced in a very wide range of environments, from circumstel- 
lar outflows to protostellar nebulae. 

Nudeation: The only area of grain research where we have 

not made significant progress over the past several years is in our 
ability to predict when refractory grains will condense from the 
vapor. Theoretical treatments of the nucleation process [9,10] are 
almost exclusively based on experimental studies of volatile com- 


pounds, such as water-vapor or low-molecular-weight hydrocar- 
bons, and are marginally applicable — at best — to formation of 
refractory materials in astrophysical environments [11]. Our recent 
efforts have been devoted to expanding the experimental database 
on the nucleation of simple, refractory materials using a single, 
experimental technique and testing various theoretical treatments 
against our results. It is obvious that grains forming in expanding, 
cooling vapors will anneal to significantly different degrees de- 
pending on the temperature where they nucleate. In all cases, 
grains formed at higher temperatures will look much more crys- 
talline than will grains formed in cooler environs [12]. ISO obser- 
vations of crystalline magnesium silicate grains in high mass loss 
rate AGB outflows have provided support for our models of both 
the growth and annealing stage of refractory evolution [13,14]. 

References: [1] Withey P. A. and Nuth J. A. (1999) Icarus , 

in press. [2] Nuth J. A. and Wilkinson G. M. (1995) Icarus , 111 ', 
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J. A. (1999) LPS XXX. [6] Rietmeijer F. J. M. and Kamer J. M. 
(1999) JCP , 110, in press. [7] Rietmeijer F. J. M. et al. (1999) Eu- 
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THERMAL EMISSION SPECTRA OF MERCURY. A E. 
Potter 1 , B. L. Cooper 2 , R. M. Killen 3 , and T. H. Morgan 4 , ^unar 
and Planetary Institute, 3600 Bay Area Boulevard, Houston TX 
77058, USA, 2 B. Cooper and Associates, Houston TX 77058, 
USA, 3 Southwest Research Institute, 6220 Culebra Road, San 
Antonio TX 78228, 4 MaiI Code SR, NASA Headquarters, Wash- 
ington DC 20546, USA. 

Introduction: The emission spectrum of Mercury may hold 

clues as to the mineral composition of the surface from reststrahlen 
bands. Unfortunately, these bands are very weak in finely pow- 
dered materials, such as are expected on the mercurian surface, 
and they lie in a region of the spectrum where there are substan- 
tial atmospheric absorptions. We have been measuring the 8- 
14-jam spectra of Mercury and the Sun at a resolution of about 
0.03 wavenumbers, using the Fourier Transform Spectrometer 
(FTS) located at the National Solar Observatory, Kitt Peak, Ari- 
zona. The objective of using such a high resolution for the mea- 
surement is to be able to look between atmospheric absorption 
lines as much as possible. The atmospheric effects are largely 
canceled by ratioing Mercury spectra to Sun spectra taken at 
nearly the same airmass, although small differences in airmass 
produce some residual atmospheric absorptions in the ratio spec- 
tra. We expect that spectral features of Mercury, if any. should 
appear in the ratio spectra, since the 8-14 pm spectrum of the Sun 
approximates a black body. 

Measurements: We have measured Fourier transform spec- 

tra of the Mercury and the Sun at intervals when we could gam 
use of the FTS at the McMath Solar Telescope at Kitt Peak. Dates 
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include January 18, 1995, April 1, 1997. November 12, 1997, and 
November 10. 1998. with a spectral range of -750-1250 wave- 
numbers. The whole planet was observed for these measurements. 
Fourier transform spectra have some unique features that must be 
taken into account when interpreting the spectra. The sampling 
frequency of the interferogram introduces a low-amplitude, high- 
frequency signal into the spectra. Slight differences in the sampling 
frequency can occur from run to run as a result of gain changes 
or slow drift in the circuits. When we divide one spectrum by 
another, these oscillations can come into exact phase with one 
another at a particular wavenumber region, resulting in a spurious 
spectral feature. We identified this effect by ratioing spectra of the 
Sun taken under different gains. Another problem has to do with 
the atmospheric absorption lines in the spectra. Where the lines are 
extremely dense, as in the ozone absorption region, and in the 
water absorption region near 8 pm, the resolution is not sufficient 
to see between the lines. Then, small differences in air mass be- 
tween the Mercury and Sun spectra can produce sufficiently dif- 
ferent absorptions in these regions to distort the ratio spectra. We 
plan to implement atmospheric corrections to these regions to 
minimize this effect. The ratio spectra are still at a resolution of 
0.03 wavenumbers, but the reststrahlen features of interest have 
widths of tens of wavenumbers. It is sensible to reduce the reso- 
lution to a value corresponding to the features of interest. Initially, 
we used Fourier and wavelet transforms to accomplish this filter- 
ing, with some success. But recently we have employed a median 
filter, which is more robust, in that it ignores wild spikes in the 
data, which the other filters tend to include. 

Results: A preliminary result from our most recent analysis 

is shown in Fig. 1. This is a measurement taken on April 1, 1997, 
when the Mercury phase was 0.562. The ratio spectrum has been 
filtered with a 5-wavenumber wide median filter. There are sev- 
eral features of interest in the spectrum. The feature near 1050 
wavenumbers results from incomplete cancellation of the ozone 
band. The feature near 940 wavenumbers is apparently an instru- 
mental effect, resulting from slight differences in sampling fre- 



Fig. 1. 


quency between Mercury and the Sun interferograms. There is a 
suggestion of the Christiansen peak beginning at about 1220 
wavenumbers. The wavenumber range is insufficient to determine 
the peak location. Finally, there appears to be a weak feature that 
might be a reststrahlen band at about 830 wavenumbers. The 
Christiansen peak and a 830-wavenumber feature have been ob- 
served by Sprague et al. [1], who suggested that a basaltic com- 
position could account for the latter feature. 

Conclusions: Mercury-Sun ratio spectra give promise for 

study of reststrahlen bands from Mercury’s surface in the spectral 
range 750-1250 wavenumbers. However, for unambiguous inter- 
pretation of these spectra, very high signal-to-noise ratios are re- 
quired to define spectral bands, with the exception of the 
Christiansen peak near 8 pm. Salisbury et al. [2] suggested that 
S/N values of the order of 1000 are required to detect reststrahlen 
bands in lunar soils. Further improvements in our data analysis 
procedures can improve our S/N past 100, but probably not to 
1000. The Christiansen peak remains the most clearly identifiable 
mineral ogical feature. 

References: [1] Sprague A. L. et al. (1994) Icarus, 109 , 156. 

[2] Salisbury J. W. et al. (1997) Icarus , 130 , 125. 
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EXES: AN ECHELON-CROSS-ECHELLE SPECTRO- 
GRAPH FOR SOFIA. M. J. Richter, J. H. Lacy. D. T. Jaffe, 
and M. K. Hemenway, Astronomy Department, University of 
Texas at Austin, Austin TX. USA. 

Introduction: We are building a mid-infrared spectrograph as 

a Pi-class instrument for the Stratospheric Observatory for Infra- 
red Astronomy (SOFIA) that will give a new perspective on how 
molecular gas turns into stars and planets. Our instrument, EXES, 
uses an echelon grating cross-dispersed by an echelle grating. 
EXES will operate at wavelengths from 5.5 to 28.5 pm and have 
a resolving power of R > 10 5 at 10 pm. At altitudes reached by 
SOFIA, the reduced water burden and atmospheric pressure will 
allow observations of H 2 , H 2 0, and CH 4 transitions that are totally 
blocked or heavily attenuated from the ground. 

Although our goals focus on high-resolution observations, we 
have designed EXES to be a versatile instrument. There will be 
three spectral modes: high spectral resolution in a cross-dispersed 
format, and medium and low spectral resolution in a long-slit 
format. There will also be a camera mode suitable for target ac- 
quisition. The low- and mid-resolution capability of EXES will 
find a broad range of applications in studies of planets, late-type 
stars, ionized nebulae, and external galaxies. 

Science Goals: We can use high spectral and spatial resolu- 
tion observations of H 2 , H 2 0, and CH 4 to study the temperature, 
density, dynamics, and chemistry of star-forming gas. H 2 rotation 
lines provide a direct tracer of warm (100 K < T < 500 K), mo- 
lecular gas such as may be found in disks, the edges of molecu- 
lar clouds, and behind shocks. H 2 0 and CH 4 will be seen in 
absorption against bright continuum sources located behind or 
embedded within molecular clouds. Doppler shifts of 30 km s _1 
are sufficient to separate the astronomical H 2 0 from telluric ab- 
sorption at SOFIA altitudes. While millimeter and submillimeter 
emission line studies of molecular clouds are biased toward the 
regions within the beam producing the strongest emission in a 
given transition, absorption line studies provide clear advantages: 
The line of sight is well-defined by the diameter of the continuum 
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TABLE 1. Typical EXES operating modes and estimated sensitivities. 


X 

(pm) 

R 

(A/ax) 

Slit Width 
(7pixels) 

Slit Length 

C) 

Bandwidth 
(pm/km s _1 ) 

NEFD 

(Jy) 

NELB 

(W rrr : sr _1 ) 

6 

120.000 

1. 4/3. 4 

7 

0.023/1200 

1.7 

3.1 x 10- 7 


14,000 

2.0/4. 8 

90 

0.023/1200 

0.6 

4.6 x 10- 7 


1800 

2.0/4. 8 

90 

0.17/9000 

0.2 

1.1 x 10- 6 

13.5 

120,000 

1. 4/3.4 

14 

0.053/1200 

4.1 

3.2 x 10-’ 


14,000 

2. 0/4. 8 

90 

0.053/1200 

1.6 

5.3 x 10- 7 


1400 

2.0/4. 8 

90 

0.53/12,000 

0.6 

2.1 x 10' 6 

18 

80,000 

2. 0/4. 8 

20 

0.13/2100 

4.6 

2.2 x 10- 7 


10,000 

2. 8/6.8 

180 

0.13/2100 

1.4 

3.6 x 10- 7 


1300 

2.8/6. 8 

180 

0.95/16,000 

0.5 

9.6 x 10~ 7 

27 

60,000 

2. 8/6. 8 

28 

0.19/2100 

7.6 

1.7 x 10- 7 


6800 

4.0/9. 6 

180 

0.19/2100 

3.3 

2.8 x 10- 7 


700 

4. 0/9. 6 

180 

1.9/21,000 

1.1 

9.3 x 10- 7 


source, so all measurements refer to the same volume of gas, and 
it is possible to observe spectral lines arising from many different 
energy levels, often at the same time. By matching our resolution 
to the width of the lines, we maximize our sensitivity to weak 
features and simplify subsequent quantitative analysis of the physi- 
cal conditions and the kinematics in star-forming gas. 

Although our own interest in the properties of molecular gas 
in star forming and protostellar regions relies almost entirely on 
EXES’s high-resolution capability, the instrument design does not 
compromise the capabilities of the low- and moderate-resolution 
long-slit mode. EXES’s resolving power ranges of -1500, 10 4 , 
and 10 5 , its high sensitivity, and its ability to observe at many 
wavelengths inaccessible from the ground make EXES on SOFIA 
a valuable tool for many types of studies including planetary at- 
mospheres, recombination lines, fine-structure lines, extragalactic 
PDRs, circumstellar envelopes, novae, and supemovae. 

Optical Design: The EXES optics consist of three sections: 

focal reducing fore-optics, an echelon chamber, and a grating 
chamber with the cross-dispersion/long-slit gratings. The fore- 
optics image the telescope secondary onto a cold pupil and 
reimage the telescope focal plane with 2:1 demagnification through 
a filter wheel onto a slit wheel. The echelon chamber contains the 
high-resolution echelon grating and its collimator/camera mirror. 
(An echelon is a coarsely -ruled grating used at very high incidence 
angle.) The grating chamber contains an echelle and a low-order 
grating, their collimator/camera mirror, focal reduction lenses, and 
the detector array. 

EXES can be used in four different modes. In high-resolution 
cross-dispersed mode (R - 10$), all three chambers are used, with 
the echelle grating serving as the cross-disperser for the echelon. 
In medium-resolution long-slit mode (R = 10 4 ), the echelon cham- 
ber is not used and the echelle serves as the primary disperser. In 
low-resolution long-slit mode (R = 1500), the low-order grating 
mounted back-to-back with the echelle is used. In acquisition 
camera mode, the low-order grating is turned face-on to act as a 
low -efficiency (r\ ~ 10%) mirror. 

Capabilities: The calculated sensitivity of EXES in the three 

spectroscopic modes, and at four representative wavelengths, is 
given in Table 1. The wavelengths selected are at typical grating 


angles. The slit widths were selected to give particular resolving 
powers. In calculating the resolving power, we included diffrac- 
tion and aberrations as well as the slit width. Of course, the echelle 
and low-order grating resolutions vary with grating angle. High- 
resolution slit lengths are chosen to avoid echelon order overlap. 
Long slit lengths are matched to the unvignetted field. Our esti- 
mate for Noise Equivalent Flux Density (NEFD, appropriated for 
observations of a pointlike continuum source) is for 10o:100s 
including nodding along slit. The estimate for Noise Equivalent 
Line Brightness (NELB, appropriate for observations of an angu- 
larly extended emission line source) is for 10a: 100s including 
nodding off slit. 

At those wavelengths where our groundbased instrument can 
observe on Gemini, the NEFD will be -5x better on Gemini and 
the NELB ~2x worse. Although the NELB is more often the 
parameter of interest for observations made by a high-resolution 
spectrometer, the gain in sensitivity that results from the lower 
telescope temperature of SOFIA is too small to justify observations 
at wavelengths with acceptable atmospheric transmission from 
Mauna Kea. Rather, we conclude that the sensitivity achievable on 
SOFIA is comparable to that achievable on the largest ground- 
based telescopes, allowing us to observe the same sources as 
would be observed by Gemini at wavelengths inaccessible from 
Mauna Kea. 

Acknowledgments: This work is sponsored by grant USRA 

8500-98-008. 
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INVESTIGATING THE MARTIAN ENVIRONMENT WITH 
THE MARS GLOBAL SURVEYOR THERMAL EMIS- 
SION SPECTROMETER. T. L. Roush, Mail Stop 245-3, 
NASA Ames Research Center, Moffett Field CA 94035-1000, 
USA (troush@mail.arc.nasa.gov). 

The Thermal Emission Spectrometer (TES) onboard Mars Glo- 
bal Surveyor (MGS) is being used to investigate the surface and 
atmosphere of Mars and the martian moons Phobos and Deimos. 
As such, it builds upon infrared observations obtained by previ- 
ous spacecraft, e.g.. Mariner 9 and Viking Orbiters [1-4]. 
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The objectives of the TES experiment are: (I) determine and 
map surface minerals, rocks, and ices; (2) study the atmospheric 
dust composition, particle size, and spatial and temporal distribu- 
tion; (3) investigate condensate clouds, C0 2 and H 2 0, location, 
temperature, and height; (4) investigate polar cap deposits, e.g., 
growth, retreat, and energy balance; (5) measure the thermo-physi- 
cal properties of surface materials; and (6) characterize the atmos- 
pheric structure and dynamics [5]. 

The TES instrument is based upon a Michelson interferometer 
and collects data in the 1700-200 cm- 1 region (-6-50 pm) at 5 
or 10 cm -1 resolution. There are also broad-band bolometric (4.5- 
100 pm) and solar reflectance (0.3-2. 7 pm) channels. The TES 
was designed to have a noise equivalent spectral radiance of 
1.2 x 10~ 8 W- 1 cm -2 sr _1 cm- 1 , corresponding to a signal-to- noise 
ratio of 490 at 1000 cm- 1 (10 mm) for a 270 K scene [5] and 
preflight data suggest a radiometric accuracy of about 1.2 x 1(H 
W- 1 cm -2 sr _1 cm- 1 [6]. In-flight observations indicate a small 
systematic calibration offset of about 1.2 x 10 -7 W -1 cm -2 sr _1 
cm -1 is present in the TES data [6]. 

MGS achieved Mars orbital insertion September 1 1, 1997, and 
entered the initial aerobraking phase [7]. MGS should have 
reached a circular orbit by early 1998. However, structural dam- 
age to one solar panel occurred during its deployment in the cruise 
phase and ultimately required an assessment of the extent of the 
damage and a much slower aerobraking period. This has delayed 
reaching the final circular orbit until March 1999 [7], 

TES and the other MGS science instruments began operating 
just after orbital insertion and continued to obtain data until No- 
vember 1998 when the TES was turned off to reduce power con- 
sumption on the spacecraft. TES data obtained during orbits 2 to 
53 have previously been published [6] and clearly illustrate the va- 
riety of scientific questions that can be addressed with them. There 
are clear spectral distinctions due to atmospheric dust, atmospheric 
condensates, and surface materials [6]; an example is provided in 
Fig. 1. 

The Mariner 9 spacecraft carried an infrared interferometer 
spectrometer (IRIS) that clearly defined a broad atmospheric dust 
absorption band between about 1300 and 200 cm- 1 [2]. The depth 
of this feature increases with increasing emission angle, indicating 



Fig. 1 . Emissivity spectra of atmospheric components and surface region 
of Mars. 



Fig. 2. Number density of the ice cloud (left), radiance (center), and 
temperature profile retrived from limb observation (left), all as a function 
of height. 


it is due to materials suspended in the martian atmosphere. The 
IRIS data have been interpreted to suggest that either a clay com- 
ponent or palagonite, a weathering product of basaltic glass, was 
responsible for the atmospheric dust signature (8,9). This absorp- 
tion is also readily observed by the TES. To date, there has been 
no mineralogical interpretation of the dust spectral signature based 
upon the TES data. However, the appearance of the emissivity 
maximum, commonly referred to as the Christiansen frequency, at 
wavelengths <8 pm, suggest that ultramafic materials cannot be 
responsible for the atmospheric dust signature (e.g., 10,11]. 

Atmospheric temperature profiles were obtained using the 
15-pm C0 2 absorption band from nadir observations. Using these 
profiles, an estimate of the surface temperature from a nearly trans- 
parent spectral region (-7.7 pm), a unit surface emissivity, and 
uniform vertical mixing, the dust optical depth can be determined. 
Dust optical depths, scaled to a 6.1-mbar surface, range from 
-0.15-0.25 during the early orbits [6]. During the onset of a dust 
storm (-orbit 51), centered near Nochias Terra, these optical 
depths increase by a factor of 2 to 4 and the spead and decay of 
the dust storm was readily monitored [6]. This clearly illustrates 
the capability of the TES data to monitor the spatial and tempo- 
ral variability of dust loading in the martian atmosphere. 

These same vertical temperature profiles can be used to monitor 
large-scale zonal atmospheric motion [6]. Based upon thousands 
of retrived profiles, evidence for a strong polar vortex centered 
near 50°N latitude is clearly indicated [6]. 

Vertical temperature profiles to higher altitudes can be derived 
from observations of the martian limbs with the TES and provide 
better vertical resolution than the nadir observations [6]. Using this 
technique, discrete water ice clouds were detected and character- 
ized [6]. An example of the indication of discrete clouds at alti- 
tudes of 45 and 53 km is provided in Fig. 2 [6]. 

The TES data have also been used to monitor the recession of 
the seasonal south polar cap [6]. The asymmetric nature of the 
polar cap recession monitored by TES is consistent with histori- 
cal observations [e.g., 12]. The radius of the seasonal cap is about 
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Fig. 3. Phobos spectrum (dashed line) compared to a mixture of Planck 
blackbody spectra (solid line). Approximate wavelengths (in micrometers) 
are given at the top. 

6° smaller, based upon the TES observations, when compared to 
previous estimates [6]. This difference may be due to the effects 
of subpixel mixing of bare ground and frost near the cap edge 
when comparing thermal and visual estimates of the cap edge [6]. 

Early interpretations of surface signatures [6] are currently 
under reevaluation due to the recognition of the significant distri- 
bution and occurrence of atmospheric condesates. Nonetheless, 
there are clear indications of spatially restricted occurrences of 
hematite on the martian surface. A region of -500 km diameter is 
centered near 1°S latitude and 3°W longitude. 

In addition to the martian surface, the TES has also observed 
the martian moon Phobos. These observations clearly indicate that 
Phobos’ surface is composed of a fine particulate powder. Initial 
interpretations suggest that the surface particulates have a very low 
thermal inertia and are likely -1 m thick. Compositional interpre- 
tation of the Phobos data is currently underway. An example spec- 
trum is shown in Fig. 3. 

References: [1] Hanel R. et al. (1972) Science , 175 , 305. 

[2] Conrath B. et al. (1973) JGR, 78, 4267. [3] Kieffer H. et al. 
(1970) Science, 193, 790. [4] Kieffer H. et al. (1977) JGR, 
82, 4249. [5] Christensen P. et al. (1992) JGR, 97, 7719. 
[6] Christensen P. et al. (1998) Science, 279, 1692. [7] Albee A. 
et al. (1998) Science, 279, 1671. [8] Toon O. et al. (1977) Icarus, 
30, 663. [9] Clancy R. et al. (1995) JGR, 100, 5251. [10] Logan 
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DEVELOPMENT OF THERMAL INFRARED EMISSION 
SPECTROSCOPY FOR GEOLOGICAL INVESTIGATIONS 
OF EARTH AND MARS. S. W. Ruff, Department of Geology, 
Arizona State University, Tempe AZ 85287-1404, USA (ruff@tes. 
la.asu.edu). 

Introduction: The Thermal Emission Spectrometer laboratoiy 

at Arizona State University was established to support the obser- 
vations of the Thermal Emission Spectrometer instrument onboard 
the ill-fated Mars Observer spacecraft. The opportunity to refly this 


instrument was made available with the advent of the Mars Glo- 
bal Surveyor spacecraft, currently in orbit around Mars. Over the 
15 years since the laboratory was established, the thermal-IR emis- 
sion technique has matured to the point that it can serve as an 
effective and reliable tool to address a wide range of geological 
questions. An overview emphasizing laboratory -based problems, 
solutions, and applications of the technique will be presented. 

Rationale for Developing Emission Technique: Infrared 

transmission, reflection, emission, and Raman techniques are avail- 
able for measuring the vibrational spectra of geologic materials in 
the laboratory. However, only the emission technique is well suited 
for performing such measurements remotely or in a field setting. 
Reflection techniques can be used to produce comparable spectra, 
but only the hemispherical reflection version strictly obeys 
Kirchhoff s law relating reflectivity and emissivity [1]. Current 
detector technology limits hemispherical reflection measurements 
to a narrower spectral range compared with emission measure- 
ments. The desire to match the spectral range offered by the TES 
instrument and to simulate the type of measurement it would make 
motivated the development of the emission technique. 

Experimental Challenges and Solutions: It has long been 

recognized that emissivity spectra are more challenging to measure 
than reflectivity or transmissivity spectra. Factors unique to the 
emission measurement include instrument self-emission, reflected 
environmental energy, temperature stability of the sample and 
instrument, and sample temperature determination. A detailed 
analysis of each of these factors has been performed and the means 
of dealing with them have been described in [2]. The result is a 
method for generating emissivity spectra of minerals and rocks that 
is efficient, accurate, and reproducible. The apparatus developed 
for this purpose is relatively simple in its design and could be 
incorporated into most commercial FT1R spectrometers. 

Spectral Library Development: Having overcome the chal- 

lenges of the emission measurement, it has been possible to de- 
velop a library of emissivity spectra. The emphasis of this activity 
has been to gather examples of common rock-forming minerals 
that are well characterized and of great purity [3]. Most of the 
samples in this collection are of a single particle size range (710- 
1000 pm) to ensure uniform spectral behavior, i.e., to avoid par- 
ticle size effects. Future additions to the library will include a 
greater range of particle sizes and the inclusion of rock spectra. 
The full library will be made available in digital form to all re- 
searchers. 

Spectral Deconvolution Capability: It has been well estab- 

lished that thermal-IR spectral mixtures can be deconvolved lin- 
early [e.g., 4]. This allows for the development of linear decon- 
volution strategies to interpret mineral mixtures. The ASU group 
has developed [4] and is actively using a linear least-squares al- 
gorithm that incorporates subsets of the spectral library to decon- 
volve rock and soil spectra into endmember components [e.g., 
5-7]. This same algorithm will be employed for the deconvolution 
of Mars surface spectra. 

Application to Mars: The Mars Global Surveyor Thermal 

Emission Spectrometer can be described simply as a FTIR spec- 
trometer with a telescope placed in orbit around Mars. It covers 
a spectral range of -1670—200 cm -1 with 143 spectral channels 
and a 3 x 2 detector array. The field of view of the detector ar- 
ray combined with the motion of the spacecraft in its polar orbit 
will produce an image cube of the entire planet at 3-km spatial 
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resolution over the 2-yr lifetime of the mission [8]. The spectral 
features of the martian atmosphere dominate the TES spectra of 
Mars, requiring atmospheric correction of the data before surface 
mineralogy can be addressed- To date, significant progress has 
been made in this effort, allowing for geologic observations of the 
planet. Some of these results will be presented during the work- 
shop. 

References: [1] Salisbury J. W. et al. (1994) JGR, 99, 
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EXTRAZODIACAL LIGHT DETECTION PROSPECTS 
WITH THE KECK INTERFEROMETER. E. Serabyn, Mail 
Stop 171-113, Jet Propulsion Laboratory, 4800 Oak Grove Drive, 
Pasadena CA 91109, USA (Gene.Serabyn@jpl.nasa.gov). 


Introduction: The Keck interferometer currently under devel- 

opment will include a midinfrared nulling experiment that is de- 
signed to search for extrazodiacal dust around nearby stars. The 
potential of this technique, as well as the basic experimental ap- 
proach, will be discussed. The initial goal is to search for 
extrazodiacal emission around nearby stars at the 10 solar-zodi- 
equivalent level, with a long-term goal of 1 solar zodi. 
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SILICON CARBIDE: THE PROBLEM WITH LABORA- 
TORY SPECTRA. A. K. Speck 1 , A. M. Hofmeister 2 , and M. J. 
Barlow 1 , department of Physics and Astronomy, University 
College London, Gower Street, London WC1E 6BT, UK (aks@ 
star.ucl.ac.uk), department of Earth and Planetary Science, 
Washington University, St. Louis MO 63130, USA. 


Introduction: Most of the solid material in the solar system 

is believed to have originated as small particles that condensed in 
outflows from evolved stars. However, most solar system solids 
(predominantly silicates) have been reprocessed and/or homog- 
enized so extensively that even the most primitive meteorite sili- 
cate samples have lost the isotopic indicators of their presolar 
origins. But some types of dust particles in the solar system have 
not been reprocessed and can potentially be associated with their 
stellar origins. One such dust type, silicon carbide (SiC), has been 
found in primitive meteorites [1] and is believed to be a signifi- 
cant constituent of the dust around C-rich asymptotic giant branch 
(AGB) stars [23], Silicon carbide grains can be divided into two 
basic groups: a-SiC if the structure is one of the many hexago- 
nal or rhombohedral polytypes, and (3-SiC if the structure is cu- 
bic [4], Silicon carbide grains exhibit a strong midinfrared feature 
between 10 and 12 pm with the peak of the p-SiC feature occur- 
ring about 0.4 pm shortward of that for a-SiC. Until now the 
observed peak wavelengths of the SiC feature in astronomical 
spectra have been interpreted as indicating a-SiC to be the domi- 
nant type of SiC around C stars [e.g., 5-9]. In fact, [8,9] found 
no evidence of (3-SiC in these circumstellar environments. Silicon 


carbide grains found in meteorites have isotopic compositions that 
imply that most of these grains were formed around C stars with 
small amounts forming around novae and supemovae [see 10-12 
and references therein]. With the exception of a few grains found 
in and among the meteoritic diamonds (and of unknown origin 
[13]) studies of meteoritic SiC grains to date have found them to 
be exclusively of the p-type [14]. P-SiC will transform into a-SiC 
above 2100°C, but the reverse process is thermodynamically un- 
likely. There is therefore an apparent discrepancy between the 
meteoritic and astronomical SiC-types, which has been discussed 
in detail by [8,9]. 

This discrepancy is addressed by discussing the problems as- 
sociated with previous interpretations of astronomical spectra that 
made use of laboratory infrared studies of terrestrial SiC. We have 
found that the problem in reconciling astronomical observations 
with meteoritic laboratory work has arisen from a misunderstand- 
ing and misinterpretation of laboratory spectra of SiC (and poten- 
tially other dust species). The results of fitting techniques used 
previously [5-9] are shown to be invalid due to the application of 
inappropriate matrix corrections to measured SiC optical proper- 
ties. Furthermore, we show that IR spectra taken using the dia- 
mond anvil cell (DAC) thin film method provides data applicable 
to astronomical situations without further manipulation. 

Laboratory Techniques and Results: New infrared (IR) 

absorption measurements of thin films of a- and (3- SiC created by 
compression in a diamond anvil cell will be presented. Unlike 
some other methods, a dispersive medium (such as potassium bro- 
mide, KBr) is not used. This relatively new approach is quantita- 
tive if sufficient care is taken to produce an appropriately thin and 
uniform film as shown by comparison of thin film spectra of vari- 
ous minerals to reflectivity data from the same samples [15,16 
and references therein]. Moreover, thin film spectra of garnets are 
nearly identical to single-crystal absorption data acquired in a 
vacuum [15], hence thin film spectra can be applied to astronomi- 
cal data without further manipulation. Our measurements strongly 
suggest through comparison of the new thin film data with previ- 
ous IR spectra collected for fine-grained KBr dispersions (in which 
the dust particles are dispersed in a KBr pellet) that the “matrix 
correction” wavelength shift invoked by [17] and adopted by other 
authors [e.g., 18,19] should not be applied to laboratory spectra of 
submicrometer grain size dispersions of SiC (see Fig. 1). 

Results of Fitting Technique: Having established that pre- 

vious fits of laboratory spectra for SiC to astronomical spectra 
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Fig. 1. Thin film spectra (solid lines) and the unconected KBr dispersion 
data [19] (dashed tines). 
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TABLE 1. Results of the X 2 -fitting of C star spectra. 


Source 

T 

1 col 

SiC Type 

T 

1 BB 

T sic 

^SiC 

T 9 7 

X'r 

IRAS 21489 + 530P 

455 

P-SiC 

449 

293 

— 

— 

0.515 

IRC + 10216* 

520 

P-SiC 

511 

230 

— 

— 

1.260 

AFGL 5076* 

525 

P-SiC 

557 

298 

0.137 

~ 

0.369 

AFGL 2494* 

525 

P-SiC 

516 

383 

0.167 

— 

0.306 

AFGL 3099* 

600 

P-SiC 

726 

329 

0.242 

— 

1.370 

AFGL 5102* 

600 

P-SiC 

650 

355 

0.161 

— 

0.345 

AFGL 2155* 

615 

P-SiC 

734 

288 

0.235 

— 

0.418 

IRAS 02152 + 2822* 

675 

P-SiC 

548 

519 

0.223 

— 

0.504 

IRC + 40540* 

680 

P-SiC 

859 

313 

0.173 


0.508 

AFGL 2368* 

800 

P-SiC 

727 

321 

— 

— 

1.772 

V Hya* 

865 

P-SiC 

1129 

393 

0.211 

— 

0.761 

IRC + 00365* 

975 

P-SiC 

1788 

215 

0.114 

— 

2.316 

CIT 6* 

1100 

P-SiC 

960 

363 

0.217 


1.294 

IRC + 50096 t 

1200 

P-SiC 

940 

455 

— 

— 

2.166 

R For* 

1400 

P-SiC 

906 

800 

— 

— 

1.237 

R Lep* 

1500 

P-SiC 

1284 

573 

— 

— 

5.267 

UU Aur* 

1500 

P-SiC 

2505 

446 

0.165 

— 

1.105 

V Cyg* 

1500 

P-SiC 

2556 

568 

0.139 

— 

1.014 

CS 776 + 

1600 

P-SiC 

993 

576 

— 

— 

1.879 

V414 Per* 

1600 

P-SiC 

1102 

920 

0.177 

— 

0.579 

AFGL 3068* 

377 

P-SiC 

394 

62 

0.030 

— 

0 092 

IARS 02408 + 5458* 

320 

P-SiC 

388 

96 

0.152 

— 

1.686 

AFGL 2477« 

340 

P-SiC 

377 

114 

0.073 

0.104 

0419 

AFGL 5625" 

333 

P-SiC 

358 

185 

0.097 

0.113 

0.306 


t 97 is the optical depth at 9.7 pm. 

* Fits with self-absorbed emission. 

1 Fits with pure emission only. 

t Fits with self-absorbed absorption. 

* Fit with Trapezium interstellar silicate, 
i Fit with p Cep interstellar silicate. 


have been erroneous due to the unnecessary application of a KBr 
correction factor, we have refitted our own UKIRT CGS3 spectra 
of carbon stars [8,9] without such a correction. We have used the 
same y} minimization routine [8], but the data to which the KBr 
correction had previously been applied [8,9] were used uncor- 
rected this time. The results in Table 1 show that there is an 
obvious predominance of the (3-SiC phase and that there is now- 
no evidence for the a-SiC phase at all. 

Conclusions: There is now no astronomical evidence for the 

presence of a-SiC in the circumstellar regions around C stars. 
While a-SiC might exist in small quantities, all observations to 
date are consistent with the exclusive presence of (3-SiC grains. 
This resolves the past discrepancy reconciling astronomical obser- 
vations and meteoritic samples of SiC grains. Furthermore, we 
have shown that laboratory IR spectra applicable to astronomical 
environments can be produced using the DAC thin film method. 
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THERMAL EMISSION SPECTROSCOPY OF SOLAR 
SYSTEM REGOLITHS. A. L. Sprague, Lunar and Planetary 
Laboratory, University of Arizona, 1629 East University Boule- 
vard, Tucson AZ 85721-0092, USA. 

Many years of difficult Earth-based observations analysis labo- 
ratory work and modeling have laid a firm foundation for what can 
be learned from thermal emission spectroscopy of regoliths in the 
solar system. At first, most observations were of rather low resolv- 
ing power and most groundbased successes lie in measuring ab- 
solute fluxes of the Moon, Mercury, the Galilean satellites, and 
asteroids through terrestrial atmospheric windows. Measurements 
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in the 1960s determined that these remote bodies have thermal 
flux characteristics indicative of silicate surfaced bodies. System- 
atic thermal spectral photometry of the Moon and many diverse 
asteroids allowed the study of the emissivity relation with phase 
and local illumination angle. More recently, spectrometers utilize 
technically advanced gratings electronics and detectors, permitting 
high signal-to-noise ratio thermal emission spectroscopy using 
resolving powers suitable for examining spectral features that are 
diagnostic of the actual chemical makeup of rocks and minerals 
comprising the regolith. 

To date, actual detailed compositional measurements of distant 
regoliths by thermal emission spectroscopic techniques are few 
because astronomers are inhibited by (1) the Earth’s atmosphere — 
it restricts the use of spectral regions to those where atmospheric 
absorptions and emissions are slight and can be removed; (2) the 
lack of availability of reference laboratory emission spectra of 
rocks and minerals (ideally obtained in simulated conditions of 
insolation, low or no atmospheric pressure, and a variety of grain 
sizes); and (3) the scarcity of facility instruments designed for re- 
motely sensing regoliths at infrared telescopes, e.g., appropriate ap- 
ertures and long-slit (near the diffraction limit and with array 
detectors for good spatial resolution), large dynamic range (to ac- 
commodate both standard stars and warm planets), suitable resolv- 
ing power (idealized for spectral activity seen in particulate mineral 
and rock mixtures). 

Actual identification of rock types on solar system bodies have 
been made by measuring the wavelength of the emissivity maxi- 
mum that is associated with the Christiansen frequency (for sili- 
cates, normally between 7.5 and 9.5 pm). Such identifications have 
been made for several locations on the Moon and a few locations 
on Mercury. Emission spectra from asteroids, the Galilean satel- 
lites, the Moon, and Mercury exhibit spectral features that are not 
yet fully understood. Actual mineral identifications resulting from 
thermal emission spectroscopy of regoliths are (1) Moon: a spec- 
tral minimum at 10.5 pm identifying olivine, and (2) Mercury: 
spectra between 7.7 and 13 pm identifying plagioclase feldspar 
with the presence of Na-rich albite at about 20-40%. 

In the future the deficiencies listed above may be overcome, 
and groundbased observing and modeling may add more successes 
to the list. In addition, the Stratospheric Observatory for Infrared 
Astronomy (SOFIA) and the Space Infrared Telescope Facility 
(SIRTF) will offer spectroscopic platforms above most (SOFIA) 
and all (SIRTF) the molecular absorptions of the Earth’s atmos- 
phere. Both DLR (Institute for Planetary Exploration, Berlin) and 
NASA (Arizona State University) have funded measurements of 
emission spectra of rock mineral and particulate samples in one 
atmosphere environments. 

■? o 

SHADOW-HIDING EFFECT IN REGOLITHLIKE MEDIA. 

D. G. Stankevich and Yu. G. Shkuratov, Astronomical Observa- 
tory of Kharkov State University, Kharkov 310022, Ukraine. 

To estimate emissivity /reflectivity properties of the planetary 
regoliths it is necessary to calculate the shadow-hiding effect for 
this class of media. Usually analytical models of the shadow-hid- 
ing effect are very cumbersome and direct computer simulations 
are more preferable. Our computer model of the effect supposes 
that (1) a regolithlike surface consists of a great number of ran- 


domly packed spherical particles; (2) surface elements of these 
particles emits as the absolutely black body (or scatters radiation 
according to the Lambert law); (3) the geometrical optics approach 
is valid; (4) the first order of scatter/emitting dominates. Modifi- 
cations of our computer model are possible, e.g., these can be 
(1) a monolayer of particles on an even substrate; (2) polydisperse 
media; (3) layered media with different statistical characteristics of 
layers. Our calculations show (1) for statistically homogeneous 
media the only significant characteristics of angular emissivity/ 
reflectivity distributions is the volume packing density (the smaller 
density the more anisotropic emission/scatter indicatrix of regolith- 
like media); and (2) maximal energy flux is observed at the zero 
phase angle for scatter and at the surface normal for emittance: for 
layered surfaces, even a thin upper layer changes the angular dis- 
tribution significantly. 

THERMAL EMISSION SPECTROSCOPY WITH THE 
SPACE INFRARED TELESCOPE FACILITY. J Van Cleve, 
Space Infrared Telescope Facility, Infrared Spectrograph Boulder 
Field Office, FA-3, Room 16321W, Ball Aerospace, 1600 Com- 
merce Street, Boulder CO 80306, USA. 

The Space Infrared Telescope Facility (SIRTF) will have an 
Infrared Spectrograph (IRS) as one of its three focal plane instru- 
ments. The IRS is composed of four modules, two of which are 
relevant to the study of dust, the solid component of disks, and 
regoliths. The Short Lo/Peak Up (SUPU) module performs spec- 
troscopy between 5.5 and 14 pm and imaging in broad bands 
centered at 15 (blue) and 22 (red) pm. Both blue and red peak- 
up fields are T x 1.3'. The Long Lo module performs spectroscopy 
between 14 and 40 pm. 

The instrument performance requirements for the IRS were 
derived in part from the science goal of establishing the miner- 
alogy, particle size distribution, and thermal profile of dusty cir- 
cumstellar disks. For example, the spectral resolution of both low- 
resolution modules is between 60 and 120, which is well-matched 
to the structure of solid-state emission features in micrometer-sized 
olivine grains. The wavelength coverage is well suited to dust at 
temperatures between 600 K and 60 K, corresponding roughly to 
the region of our solar system between Venus and Uranus. 

The same instrument capabilites can be used to study absorp- 
tion and emission features in interstellar dust in our own galaxy 
and identify the composition and temperature of dust in other 
galaxies. The redshift of distant galaxies can be determined from 
these features as well. Solar system science enabled by the instru- 
ment includes minerology of cometary grains, identification of 
solid hydrocarbons such as PAHs or tholins, and measurement of 
the thermophysical properties of main belt asteroids, satellites of 
the outer planets, and loose icy bodies between the orbits of Ju- 
piter and Uranus. For objects colder than 60 K, SIRTF observa- 
tions with the spectral energy distribution (SED) of the Multiband 
Imaging Photometer (MIPS) would be useful for thermal studies. 

DUST EMISSION FROM HERBIG Ae/Be STARS: EVI- 
DENCE FOR DISKS AND ENVELOPES. D Vinkovic 1 , A 
Miroshnichenko 2 , Z. Ivezic 3 , and M. Elitzur 1 , department of 
Physics and Astronomy, University of Kentucky, Lexington KY 
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40506, USA (dejan@pa.uky.edu; moshe@pa.uky.edu), 2 Pulkovo 
Observatory, St. Petersburg 196140, Russia (amiroshn@infopro. 
spb.su), department of Astrophysical Sciences, Princeton Uni- 
versity, Princeton NJ 08544, USA (ivezic@astro.princeton.edu). 

Introduction: To date several attempts have been made to 

determine the geometry of dust around Herbig Ae/Be (Haebe) stars 
from IR and millimeter-wave data. The analysis was based on one 
of two possibilities: an optically thick disk, or an optically thin 
spherical envelope. From the shape of the spectral energy distri- 
bution (SED), Hillenbrand et al. [1] suggested that IR emission 
from 30 Haebes is dominated by optically thick accretion disks, 
while Hartmann et al. [2] proposed envelopes. Both models repro- 
duce the observed XF x <^k- m shape. 

Miroshnichenko et al. [3] were able to successfully fit the 
detailed SED of eight Haebe stars with optically thin evelopes in 
spherical free fall. The motivation was that the free fall model, like 
a Keplerian disk, gives XFx~X- 4/3 . With this model we can explain 
the strong silicate 10-pm feature in emission, which indicates 
optically thin dust. However, imaging at millimeter- wavelengths by 
Mannings and Sargent [4] (MS hereafter) clearly shows disks 
around their sample of seven Haebe stars. They also find that thin 
envelopes cannot explain the observed image sizes. 



Fig. 1. Modeling of the Haebe stars in MS sample. Short-dashed line is 
the SED of star and envelope, long-dashed line is the disk’s SED, and 
solid line is the total theoretical SED. 


Furthermore, Henning et al. [5] were able to fit the SED of AB 
Aur by spherical envelope, while its observed small size rules out 
such dust configuration. They use this contradiction to rule out 
spherical models even though AB Aur shows the silicate 10-pm 
feature in emission. Additional controversy is introduced with 
imaging at far IR wavelengths by Di Francesco et al. [6]. They 
have observed objects with image size at 100 pm smaller than the 
one at 50 pm. This contradicts the fact that the temperature de- 
creases with distance from the central star and therefore the im- 
age size at longer wavelengths must be larger. 

We propose a simple solution that resolves the apparent incon- 
sistencies of observations at different wavelengths, namely, the 
dust distribution in these sources has both spherical and disk com- 
ponents. The former can be large and dominates the SED at IR 
wavelengths; the latter is smaller and dominates at longer wave- 
lengths. Because of its very low surface intensity, the envelope is 
not visible at millimeter-wavelengths but becomes detectable and 
bright at shorter wavelengths. 

The Model: The fact that there are indications of optically 

thin and thick dust in the same object leads us to a model where 
a thick disk is surrounded and heated by a thin envelope. Such 
heating is possible because the envelope has higher temperature 
than the disk at the same distance from the star. The reason for 
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Fig. 2. Modeling of the Haebe scars in MS sample (see Fig. 1 for details). 
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having a cooler disk is twofold: It is heated by the central star at 
a grazing angle, and it is optically thick. 

We model the radial density profile of the envelope with a 
broken power law. It extends from an inner radius corresponding 
to dust evaporation (we use an evaporation temperature of 1500 K) 
to an outer radius R out . 

We consider also some transitional radius Rp where the power 
p of the density power law p(r)~r~P changes from p in to p ou{ . The 
dust in our models consists of the standard ISM mixture of sili- 
cate and graphite grains from Draine and Lee [7]. The disk is 
geometrically thin and has an inner radius at the stellar surface 
with temperature 2~ 1/4 T*, where T* is the stellar temperature. It ex- 
tends to a radius R^, specified by the disk’s outer temperature T d . 
The transition of the disk from optically thick to thin is defined 
by x 350 , the optical depth at 350 pm. When t 350 = 1, the disk will 
be thin at X > 350 pm. Modeling with a thick disk is plausible 
when the optical depth at the wavelength corresponding to T d is 
larger than 1. 

We have analyzed this model with the aid of the scaling theory 
by Ivezic and Elitzur [8], where the radius is dimensionless Y = 
R/Rj (Rj is the inner envelope radius). The SED was calculated 
by the code DUSTY [9]. The stellar spectral shape was taken from 
Kurucz [10]. 

Results and Conclusions: We applied the described model to 

the MS objects. In our search for the fitting parameters we used 
0.5 < p in < 2, 50 < Y p < 500, 0 < p out < 2. 150 < Y ou( < 10000, 
x 350 > 0.2, and T d < 50 K. The optical depth of the envelope was 
described by the visual optical depth t v . The correction for inter- 
stellar extinction A v was considered a free-fitting parameter. It is 
mainly constrained by the short-wavelength data and results are 
consistent with the published estimates. The fitting yields a range 
of possible parameters for each object, and in all of them we are 
able to satisfy the observed size constraints and fit the SED from 
UV to millimeter-wavelengths (see Figs. 1 and 2). We can also 
rule out the flared disk model described by Chiang and Goldreich 
[11] because their optically thin layer surface corresponds to our 
envelope with a much smaller t v - 0.01 than we get from the fit 
(t v is always larger then 0.1 and smaller then 1 except in CQ Tau, 
where it is —2.5). The near IR is dominated in almost all objects 
by the envelope with 1.5 < p in < 2. We also give an upper limit 
on T d for each object. The inclination angle is not a free param- 
eter and can be derived from the fitting results. Usually it has a 
large uncertainty. The only noticeable difference between the 
models and observations is in the 10-pm region, and may be due 
to somewhat incorrect chemical composition and/or gram size dis- 
tribution. Details of our results and theory will be published in a 
separate publication. 

We conclude that the dust around MS objects includes both an 
optically thick and compact disklike component, and an optically 
thin, more extended component. The disk extends to several ICHR* 
(R. is the stellar radius). The envelope starts at -100R, and ex- 
tends to -lC^R*. Dust density in the inner envelope (up to several 
KPR* similar to the disk extent) resembles a steep power-law- 1/ 
r 2 and flattens out to roughly 1/r at larger radii. 
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THERMAL EMISSION SPECTROSCOPY OF LABORA- 
TORY REGOLITHS. C. Wagner, Institute of Planetary 
Exploration, German Aerospace Center (DLR), Rudower Chaussee 
5, Berlin, Germany. 

Introduction: Thermal emission spectra of samples of known 

mineralogic composition and well-defined grain size and porosity 
are valuable for interpreting remote sensing emissivity data ob- 
tained for the various planetary objects in the solar system. Be- 
cause spectral reflectance of a material is more convenient to 
measure in the laboratory than spectral emissivity, it is common 
practice to interpret emissivity spectra of planetary objects by use 
of reflectance spectra of terrestrial minerals or rocks. The then 
required conversion of spectral reflectance (R) at any given wave- 
length into spectral emissivity (E) at this wavelength is carried out 
by a complementary relationship between emissivity and reflec- 
tance E = 1 - R according to Kirchhoff s law. But Kichhoff s law 
holds only under conditions of thermodynamic equilibrium, and is 
violated when the emitting material is nonisothermal. Planetary 
objects with little or no atmosphere, however, can evolve signifi- 
cant thermal gradients in the uppermost surface layer due to near- 
surface radiative cooling [1,2]. The effect of thermal gradients on 
the emissivity spectra of particulate minerals is complicated, and 
unfortunately “the spectral emissive properties of particulate ma- 
terials in nonterrestrial conditions are understood neither theoreti- 
cally nor experimentally” [3]. At any rate, use of reflectance data 
can provide only a more or less justified approximation for the 
actual emissivity spectrum. The (unknown!) degree of approxima- 
tion depends on the thermal gradient within the radiating surface 
layer and is difficult to estimate for a remote target. This gives 
reasons for measurements of spectral emissivity of minerals and 
rocks in a simulated space environment (i.e., at low atmospheric 
pressure). But up to now such measurements are very scarce and 
comprise (as far as I know) only six published studies [1,3-7]. 

An important point is that even if emissivity measurements 
carried out at atmospheric pressure cannot provide a true analog 
for emissivity spectra of airless planetary objects, they still prove 
to be a valuable and interesting link and intermediate result be- 
tween (1) emissivity measurements in vacuo and (2) reflectance 
measurements at atmospheric pressure. This is due to the 
following: For a material with a (vertical) temperature gradient the 
emission spectrum is controlled not only by the materials emissiv- 
ity, but also by the interplay of (1) the wavelength-dependent 
opacity of the material, and (2) the varying brightness temperature 
[cf. 8-11]. 

Experimental Procedure: The apparatus and experimental 

procedure are similar to that used by Christensen and Harrison 
[12] and are described in more detail in [13]. For emissivity 
measurements the sample is placed in an Al cup (3 -cm diameter) 
and heated from below to a constant temperature of 90°C. The 
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samples’ thermal radiation emitted normal to the sample surface 
within a circular area of 2.5-cm diameter is collected by a Au- 
coated parabolic off-axis mirror and reflected to the entrance port 
of a Fourier transform infrared spectrometer (Bruker IFS 88). The 
emissivity of the sample is obtained by ratioing the background- 
corrected signals obtained for the sample and for the reference 
blackbody [e.g., 3]. For a thick sample of particulate material, 
however, the vertical thermal gradient is larger than for the black- 
body (a 6-mm-thick Cu disk with a Nextel Velvet coating). There- 
fore, a systematically too low emissivity would be obtained. Accord- 
ing to common practice [3,5,10,12], this systematic error is cor- 
rected by a fit procedure. It is assumed that the sample approached 
a known maximum emissivity (usually 1 .0) somewhere in the mea- 
sured spectral range. At this wavelength a Planck curve is fitted 
to determine the actual sample temperature and then the emissiv- 
ity of the sample is scaled to that of the blackbody [12]. Fortu- 
nately, for many minerals the assumption E = 1 is justified at the 
materials’ Christiansen wavelength, i.e., where n(X) = 1. Another 
approach to get the samples’ true emissivity spectrum is keep- 
ing the temperature gradient low within the sample by using a 
very small sample thickness. Then no fit procedure is required. 
The results of the two methods practically agree as can be shown 
by experiment [14]. Using a thin sample may be advantageous if 
only a small quantity of material is available or when no fitting 
point for E is known. Though in the thm-layer approach errone- 
ous spectra are possible due to an insufficient optical depth, but 
they result only if the sample thickness falls below a certain mini- 
mum value. 

Results: Emission spectra in the wavelength range from 4 to 

17 pm were measured for the well-characterized materials listed 
below for several (usually seven) grain-size fractions ranging from 
<25 pm to 180-200 pm, and, where samples are available, spec- 
tra were also measured for the massive material. 

1. Carbonates and sulphates. Calcite magnesite siderite dolo- 
mite/ankerite; anhydrite and kieserite. 

2. Mars soil analogs and ingredients. Four palagonites and the 
JSC Mars-1 martian soil simulant (these spectra are not for sev- 
eral grain sizes!), two hematites, two maghemites, three magne- 
tites. 

3. Feldspars. Six plagioclase feldspars (albite, oligoclase, 
andesine, labradorite, bytownite, anorthite), five alkali feldspars 
(albite, anorthoclase, sanidine, orthoclase, microcline), and three 
feldspathoids (leucite, nepheline, sodalite), all for at least seven 
grain sizes. 

4. Pryoxenes. Augite, bronzite, and hedenbergite, all for the 
full range of grain-sizes. 

In the future we plan to measure emissivity from: 

1. Lunar soils of different exposure age obtained from the 
Apollo 12, 14, and 16 missions, i.e., samples 12023,83, 12070,12, 
14259,602, 60051,32, and 61241,65. 

2. Feldspar-pyroxene-mixtures. (1) 25 and 50 wt% bytownite 
in hedenbergite (i.e., Ca-rich feldspar mixed with Fe-rich pyrox- 
ene as a model for the Moon); grain sizes <63 pm and 125- 
180 pm. (2) 25 and 50 wt% albite in bronzite (i.e., Na-rich 
feldspar with Mg-rich pyroxene as a model for Mercury); grain 
sizes as in (1). 
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SOLID-STATE SPECTROSCOPY OF CIRCUMSTELLAR 
DUST: INFRARED SPACE OBSERVATORY RESULTS. R 

Waters, University of Amsterdam, Kruislaan 403, Amsterdam NL 
1098 SJ, The Netherlands. 

The Short and Long Wavelength Spectrometers (SWS and 
LWS) onboard the Infrared Space Observatory (ISO) have yielded 
a wealth of new information concerning solid-state spectra of cir- 
cumstellar dust. For the first time, a full inventory of solid-state 
bands over a wide wavelength range with moderate spectral reso- 
lution has been obtained. An overview will be given of highlights 
of the ISO results in this area for both the O-rich and the C-rich 
dust chemistry. A comparison between dust surrounding young 
and evolved stars is made with some surprising analogies. 
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THERMAL EMISSION SPECTROSCOPY OF 1 CERES: 
EVIDENCE FOR OLIVINE- F C Wittebom 1 , T L Roush 2 , 
and M. Cohen 3 , ^ETI Institute, Mail Stop 245-6, NASA Ames 
Research Center, Moffett Field CA 94035-1000, USA, 2 Mail Stop 
245-3, NASA Ames Research Center, Moffett Field CA 94035- 
1000, USA, 3 Radioastronomy Laboratory, University of California, 
601 Campbell Hall, Berkeley CA 94720, USA. 

Thermal emission spectra of the largest asteroid 1 Ceres ob- 
tained from the Kuiper Airborne Observatory [12] display features 
that may provide information on its surface mineralogy. The plot 
in Fig. 1 is the Ceres spectrum (calibrated using a Boo as a stan- 
dard) divided by a standard thermal model (STM). Also shown is 
the emissivity spectrum deduced from reflectivity measurements by 
Mustard and Hays [3] for olivine grams <5 pm in diameter. The 
general shape of the Ceres and the olivine curves agree in essen- 
tial details, such as the maxima from 8 to 12 pm, the minimum 
between 12 and 14 pm, the broad peak near 17.5 pm, and the 
slope beyond 22 pm. (Use of the 10- to 15-pm grain reflectivities 
[3] provides a better match to the 12- to 14-pm dip. We used a 
value of unity for |3, the beaming factor associated with small-scale 
surface roughness in our STM. Adjustment of (3 to a lower value 
raises the long-wavelength side of the Ceres spectrum, providing 
an even better match to the olivine curve.) The emissivity behav- 
ior roughly matches the emission coefficients calculated by Mukai 
and Koike [4] for olivine particles with a particle radius of 3 pm. 
Their calculations show not only the negative slope from 23 to 
25 pm, but a continued decrease past 30 pm. The Ceres emissiv- 
ity is thus similar to that of small olivine grains from 8 to 30 pm, 
but olivine’s emissivity is lower from 5 to 8 pm. 
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Fig. I. Ceres emissivity vs. wavelength. Measured Ceres emissivity 
values (points with error bars) are compared with those measured for small 
olivine grains by Mustard and Hays [3] denoted by an x. 


Although no olivine features have been identified for Ceres 
from reflectance spectra (0.4-2. 5 pm), such features have been de- 
tected in several small asteroids [5]. The largest of the asteroids, 
Ceres is a prime candidate to have undergone differentiation. 
Olivine is expected to be a major product of differentiation and 
typically forms at depth due to its relatively high density. Subse- 
quent impacts are expected to have emplaced the olivine-rich 
material at the surface of the differentiated asteroids and perhaps 
spread it to others as well [6,7]. The reflectance spectra of Ceres 
have been interpreted to indicate a surface consisting of hydrous 
silicate(s) with low ferric iron content and an opaque phase [8a). 

Thus we are left with the question of why reflection spectros- 
copy provides no indication of olivine on Ceres, whereas the ther- 
mal observations reported here clearly suggests its presence. 

Reflectance spectra probe the uppermost layers of a surface and 
spectral features are readily subdued or masked by trace amounts 
of opaque materials [e.g., 8b]. Surfaces emit energy from a vari- 
ety of depths, and at these wavelengths the vibrational fundamen- 
tals are so strong that they are not as readily subdued or masked 
by overlying materials [8c]. One possibility that could account for 
the absence of olivine absorption features in the Ceres reflection 
spectrum is a thin layer (say a few micrometers or less) of dust en- 
riched in opaque material. In fact, a mixture of olivine with some 
opaque material might explain the additional emissivity shortward 
of 8 pm seen in the emissivity spectrum of Ceres. While small 
asteroids are thought to be depleted in small dust grains [8], Ceres 
is large enough to retain fine-grained dust. For example, Le Bertre 
and Zellner [9] conclude that Vesta’s surface has particles larger 
than 50 pm coated with particles smaller than 10 pm. Masking of 
olivine features in the Ceres reflectance spectrum could have 
broader implications because the apparent paucity of olivine-nch 
asteroids is regarded as a mystery [6,7]. This raises the ques- 
tion: Could a fine dust layer obscure the olivine component of 
large asteroids in the near infrared, thus accounting for the appar- 
ent deficit of olivine in the asteroid belt? 

References: [1] Wittebom F. C. et al. (1999) Astrophys. 

J in press. [2] Cohen M. et al. (1998) Astrophys. J., 115, 1671— 
1679. [3] Mustard J. F. and Hays J. E. (1997) Icarus , 125 , 1 45— 


163. [4] Mukai T. and Koike C. (1990) Icarus, 87, 180-187. 
[5] Cloutis E. A. (1993) IPS XXIV, 317-318. [6] Bell J. F. et al. 
(1989) in Asteroids II (R. P. Binzel et al., eds.), pp. 921-945, 
Univ. of Arizona, Tucson. [7] Burbine T. H. (1994) Meteoritics, 
29, 453. [8] Dollfus A. et al. (1989) in Asteroids II (R. P. Binzel 
et al.. eds.), pp. 594-616, Univ. of Arizona, Tucson. Ceres emis- 
sion: Wittebom Roush and Cohen [8a] Gaffey M. J. and McCord 
T. B. (1979) in Asteroids (T. Gehrels, ed.), pp. 688-723, Univ. of 
Arizona, Tucson. [8b] Clark R. N. (1983) JGR, 88, 10635-10644. 
[8c] Christensen P. R. and Harrison S. T. (1993) JGR , 98, 19819— 
19834. [9] Le Bertre T. and Zellner B. (1980) Icarus, 43, 1 72— 
180. 

ABUNDANT COOL MAGNESIUM-RICH PYROXENE 
CRYSTALS IN COMET HALE-BOPP. D H Wooden, Mail 
Stop 245-3, NASA Ames Research Center, Moffett Field CA 
94035, USA. 

Modeling of the observed dust emission from Comet Hale- 
Bopp over a large range of heliocentric distances (2.8 AU - 0.93 
AU - 1 .7 AU) led to the discovery of Mg-rich pyroxene crystals 
in the coma. These pyroxene crystals are apparent in the 10-pm 
spectrum only when the comet is close to perihelion (r h = 1 .2 AU) 
because they are cooler than the other silicate minerals. The py- 
roxene crystals are cooler than the other species because of their 
high Mg-content: They do not absorb as efficiently as the other 
silicate minerals. Given the same high Mg content of Mg/(Mg + 
Fe) = 0.9, radiative equilibrium computations show that pyroxene 
crystals are expected to be 150 K cooler than olivine crystals. The 
pyroxene crystals are also about lOx more abundant than the other 
silicate mineral species. Their high Mg content and relatively large 
abundance are in agreement with the preponderance of pyroxene 
interplanetary dust particles (IDPs) and the recent reanalysis of the 
PUMA-1 flyby of Comet Halley. Before Hale-Bopp, only olivine 
crystals were detected spectroscopically in comets, probably be- 
cause the pyroxene crystals are less optically active, hence signifi- 
cantly cooler and harder to detect in contrast to the warmer silicate 
species. Determining the relative abundances of silicate minerals 
depends on their Mg contents. 

If the pyroxene crystals in Comet Hale-Bopp are solar nebula 
condensates, then they probably had to form during the early FU 
Orionis epoch when the inner disk was hot enough and be trans- 
ported out to the region of formation of icy planetesimals without 
being reheated. Reheating events appear to reincorporate Fe back 
into the crystals or form Fe-rich rims, which are not consistent 
with the high-Mg-content crystals. The condensation of Mg-rich 
pyroxene crystals is expected from solar nebula thermal equilib- 
rium computations. However, their subsequent transport to the 
outer solar nebula unaltered has yet to be theoretically demon- 
strated. 

The discovery of Mg-rich crystals in Comet Hale-Bopp and in 
AGB stars opens the possibility that these crystals are relic inter- 
stellar grains. One-third of IDPs have been shown to have signifi- 
cant deuterium enrichments, thus indicating that they contain 
presolar material. By spectroscopic analogy to IDPs the Mg-rich 
pyroxene crystals in Comet Hale-Bopp may be presolar grains. If 
so, then the comet contains largely ISM silicates. ISM grains may 
have been the dominant source of dust in the outer early solar 
nebula. 



